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1 Introduction 

It is well known that the spectra of integrable two dimensional QFT's in finite volume can be studied 
by the thermodynamic Bethe ansatz (TBA) approach [1]. TBA results in a system of non-linear 
integral equations which can be always "miraculously" recast into the universal Y-system — a system 
of functional equations on Y-functions, looking simpler than TBA, but at the price of loss of a certain 
analyticity input. Apart from some exceptional cases (Lee- Yang model, SS-model, Wess-Zumino 
cosets) both the TBA system and the Y-system are infinite. However, Destri and de Vega noticed [2] 
that the infinite TBA system for the vacuum of the Sine-Gordon model can be reduced to only one 
non- linear integral equation, and this was later generalized to excited states [3]. For the models with 
higher rank groups of symmetries a similar reduction to a few integral equations is usually possible. 
DdV-type equations often make possible a systematic study of IR and UV limits. In general this is 
the best one can do with the spectral problem of a finite volume integrable QFT. This reduction, still 
not very well studied and understood, can be traced down to the fact [4, 5] that the Y-system itself 
is a gauge invariant version (see (2.3)) of the integrable finite difference bilinear Hirota equation, or 
T- system: 

The equation is essentially the same for all known integrable 2D QFT's^, and what differs from 
one model to another is its boundary conditions w.r.t. the "representational" variables a, s and the 

^with slight modifications for the algebras other than Aj.. 



- 1 - 



analyticity properties w.r.t. the spectral parameter. In practice, the Hirota equation is often a very 
good starting point for attacking the finite volume spectrum problem in any integrable 2D QFT's (e.g. 
for an important case of the principal chiral field [4, 6]): even the analyticity properties of T-functions 
appear to be considerably constrained by the structure of Hirota equation. 

All these observations should hold for the AdS/CFT Y-systems [7], and in particular for the 
most studied AdS5/CFT4 case, solving the problem of planar spectrum of anomalous dimensions in 
N=4 SYM theory (see for the recent review [8], and in particular the chapters [9, 10] therein). The 
AdS5/CFT4 Y- and T-systems with T-hook boundary conditions proposed in [7] and summarized 
in figure 1 were later shown to be equivalent, with certain analyticity requirements [11-13], to the 
TEA equations [14-16]. It was shown in [17, 18] (see also [19, 20]) that the T-system, and hence the 
Y-system, in T-hook can be formally solved in terms of Wronskian determinants of a finite number of Q- 
functions — a generalization of Baxter's Q-function. But the corresponding finite system of nonlinear 
integral equations (FiNLIE), a remote analogue of Destri-de Vega equations, was still missing. In this 
paper we propose such a FiNLIE for the planar AdS/CFT spectrum. 

The difficulty of finding FiNLIE — possibly the ultimate step in simplifying the spectral problem in 
AdS/CFT — resides in its formidably complicated analytic structure w.r.t. the spectral parameter u: 
the Y-functions live on a Riemann surface with infinitely many sheets connected by an infinite system 
of Zhukovsky cuts (we will call them Z-cuts) parallel to the real axis, with the fixed branch points 
at M = ±2g + |Z. This cut structure can be already seen in the asymptotic, infinite length solution 
of the Y-system [7] deduced from the asymptotic Bethe ansatz equations [21] describing the infinite 
volume spectrum. This analytic structure is also visible from the dispersion relation of elementary 
excitations and their bound states in the string sigma-model in the light-cone gauge, when the energy 
and momentum are parameterized by the spectral parameter u (see eqs.(2.7)-(2.8)). 

One way to overcome these difficulties with analyticity is to work in the TBA formulation of 
the Y-system which implicitly includes all the information about the analyticity properties of the Y- 
functions. This approach made possible the first numerical computation of the anomalous dimension 
of the lowest lying, Konishi state in planar A/" = 4 SYM [22] . 

However, to write FiNLIE we need to explicitly understand the missing analyticity properties of 
all Y-functions or T-functions. In the papers [11, 13] the basic analytic properties of Y-functions 
were decrypted from TBA. Alternatively, together with the Y-system equations, these properties were 
shown to be equivalent to the TBA equations. This was an important progress though the resulting 
formulation of the problem contained the same infinite number of Y-functions and it still remained 
quite complex. 

We find in this paper that the analytic properties become considerably simpler and more natural 
when formulated in terms of T-functions. We managed to formulate here a simple and sufficient set 
of analyticity properties of T-functions, mostly following from the symmetries of the model. In order 
to formulate these properties in a simple way we also have to make a good choice of labeling of the 
sheets of the Riemann surface. In this paper, we describe a choice which we call "magic" , additional 
to the "physical" and "mirror" labeling already extensively used in the literature. The magic sheet 
has, unlike the mirror sheet, "short" cuts between the pairs of branch points u = zt2g + |Z and 
coinciding with the mirror sheet in the vicinity of the real axis. On this sheet, we discover a new 
quantum symmetry (already mentioned in [18]), which generalizes the classical Z4 symmetry of the 
string sigma model on the so(i^)x'sO(6) ^'^^^^ [^3, 24]. Using this symmetry, we identified certain 
Q-functions having one single short Z-cut, which means that they can be written in terms of a density 
with a finite support [—2(7, 2g]. 

These observations allow us to formulate the first (probably still perfectible) version of FiNLIE 
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and even to perform its numerical study for the Konishi dimension. The results of this numerical 
study are in perfect correspondence with the results obtained in [22, 25] from TBA and thus confirm 
the validity of our FiNLIE. We also prove the equivalence between FiNLIE and TBA analytically. 




S S 



(a) lattice for T-functions' (a,s) indices (b) lattice for Y-functions' (a,s) indices 



Figure 1: Boundary conditions for Y-system (b) and T-system (a) of AdS5/CrT4: T- 

shaped "fat hook" (T-hook). We will often distinguish in the T-hook (a) the upper band - all the 
nodes with a > \s\, the right band - all the nodes with a < s and the left band - all the nodes with 
a < \s\,s < 0. 

This paper is composed as follows: section 2 introduces the Y-system and our notations. The 
readers familiar with the subject can go directly to section 3 where the analytic properties of T- 
functions are described and the Z4 symmetry of the classical string is generalized to a new symmetry 
of the full quantum T-functions. In section 4, the Wronskian solution of T-system in terms of a finite 
set of Q-functions is introduced. We work out here a compact and general representation of such 
solutions in terms of exterior forms. In section 5, the equations constraining these Q-functions are 
obtained with the use of the analyticity and symmetry constraints identified in section 3. Finally, 
section 6 summarizes our construction and the resulting FiNLIE, while section 7 presents the results 
of our first numerical implementation of this FiNLIE. A systematic study of the analyticity properties, 
as well as a recasting of the spectral problem into the FiNLIE form, are performed explicitly for the 
sl{2) sector's symmetric states, and in particular having only two magnons, containing a relatively 
large class of operators including Konishi operator. 



2 Y-system generalities 

For the detailed description of the Y-system for AdS/CFT we refer to [8-10]. In this section we list 
the main results and general facts about Y-systems and Hirota equation. The reader familiar with the 
basics of Y-systems can skip this section. 

2.1 AdS/CFT Y-system 

The AdS/CFT Y-system belongs to a class of functional equations arising in spectral problems of 
many integrable models 

(1 + l/ya+l,s("))(l + l/Ya-l,s{u)) 

In particular, the same equations are used to compute the spectrum of SU (N) principal chiral model 
[4, 6] which is an important source of inspiration for the AdS/CFT case. Y-functions are associated 
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with certain nodes on a two-dimensional (a, s) lattice. The main difference from model to model is 
the shape of the domain on this lattice where the Y-functions are defined. In the present AdS5/CFT4 
it is the T-shaped hook [7] depicted in figure lb. ^ 

In addition to the shape of the (a, s) domain one should also specify the analyticity properties of 
the Y-functions. These properties can be deduced from the TBA equations [14-16] and were explicitly 
written in [9, 11, 13]. 

Y-functions are natural objects in the TBA approach, however not the most elementary ones. 
More fundamental quantities are the T-functions defined by 

Ya, = l^'-^'l^''-' . (2.2) 

a+l,s-^ a— l.s 

As it is the case for many relativistic sigma models, T-functions can sometimes be identified with 
eigenvalues of the transfer matrices of an underlying spin chain-type discretization. A similar scenario 
may exist in the AdS/CFT case, however the physical interpretation of the AdS/CFT T-functions and 
their operatorial realization is not yet established. In the strong coupling limit the T-functions were 
identified with the psu(2,2|4) characters of the monodromy matrix [17]. This classical monodromy 
matrix was perturbatively quantized, up to two loops on the world sheet, in [29, 30] and the Hirota 
functional equation (1.1) was shown to hold at least in this approximation. 

It was noticed for the AdS/CFT system that T-functions are much simpler objects than the Y- 
functions in all limiting cases where the general solution of the Y-system was obtained explicitly, i.e. 
in the asymptotic large volume limit [7] and in the strong coupling limit [17]. It is very natural to 
assume that the formulation of analytic properties at the level of T-functions should be also more 
revealing and simple. Indeed this is the CclSG, clS described in the next sections. 

As one can see from (1.1), the T-functions live in an enlarged T-hook depicted on figure lb. The 
meaning of the shape of the (a, s) domain in this case also becomes clear in terms of the T-functions: 
the super- Young diagrams of the unitary highest weight representations of psu(2, 2)4) should be en- 
tirely contained inside this domain [17, 27]. Furthermore, as we explain in this paper, an important 
Z4 symmetry of the super-coset model noticed in [18] has a transparent realization in terms of the 
analyticity of T-functions w.r.t. both the spectral parameter u and the representational variables 
(a, s). 

It is easy to see that the definition of T-functions by (2.2) is ambiguous. This ambiguity is 
expressed in the "gauge" transformation 

Ta.,s ^ 9\'''^'^g^^''-'^9t'-'^g[-''^'Ks (2.3) 

which does not affect Y's and the Hirota equation. Here and everywhere in the text we use a short-hand 
notation for the shift of the spectral parameter: 

/[±'^l=/(u±za/2), (2.4) 

It is also useful to notice that due to the Hirota equation (1.1) there are several equivalent ways for 
the Y-functions to be expressed in terms of the T-functions 



7^+ 7^— 7^+ 7^— 

l + l/Y^,^ ^ °'V" ■ (2-5) 



a+l,s-^ a— l.s a,s-t-l-' a,s — 1 



^More general Y-systems with T-hook shapes are possible for generalizations to other non-compact supersymmetry 
groups, see [26-28]. 
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2.2 Notations 

Let us introduce some important notations used through the text. We define the Zhukovsky variable 

47r 



X + 1/x = it/(7, where g = ^ is related to the 't Hooft coupling A, choosing two different branches'^ 




The first choice x{u) will be referred to as the mirror branch, whereas x{u) will be called either 
"physical" or equivalently'* "magic". x{u) has a long cut, called from now on the Z-cut, connecting 
the branch points ±25 through the infinity, whereas x{u) has a short cut [— 2g, 2g] called Z-cut (see 
figure 2). 

The asymptotic energy of a single magnon (or of a bound state of a magnons) is given by [31-33] 



2ig 2ig 

Similarly, the asymptotic magnon momentum is 



eaiu) = a+ , , . — —. — r . (2.7) 



I X 

The energy and the momentum of a state are given by 

M 



1 

p,(w) = ^log^. (2.8) 



S=^ei(Mfe)+^ / ^.duta{u)\og{l+Ya,Q{u)) , 

fc=l a=l-^^°° 
A/ 00 „oo 

P=^pi(wfe)+^ / ^duPa{u)\0g{l + Yafi{u)) 

1 1 1 00 



(2.9) 



k=l a=l' 



2ig 2ig . x 1, a^t+'^l 

ZTT:;! - -TTTT : Pa(w) = 7log- 



The anomalous dimension is defined as 7 = — M. The shift of the spectral parameter is always 
defined by a path which avoids the fixed Z- and Z^-cuts of the functions^. For instance, the path in 
the shift w — >■ M + ia/2 of i'"! never crosses the Z-cut, while the path in the shift of x'^I never crosses 
the Z-cut. As a consequence, ^ x^~^^ whereas x^+^l = x^+^l (see figure 2) and thus the choice of 
the cuts affects the functional relations. The "hat" symbol over a function of the spectral parameter 
u (for instance in the definition (2.6) of x{u)) means that we choose a Riemann sheet with only short 
Z-cuts. When the "hat" appears over a function of u and of the (a,s) labels, as on figure 1, then it 
denotes also the analytic continuation in the labels (a, s) performed on the sheet having only Z-cuts. 
This will be introduced in more detail in section 3.4 and section 3.5 (see for instance (3.17)). 

We will also use in this paper a short-hand notation for long cuts shifted by the distance in/ 2 
from the real axis as and similarly for the short cuts as Z^. Another important notation which 
we will frequently use is An, denoting the class of functions analytic^ inside the strip |Imu| < ^. The 



^To define the branches, we always use the definition of square roots which has a cut on R~ and is positive on R"*". 

^ The function x is quite particular because the magic and physical sheet coincide. We will see that there exists other 
functions (like T- and Y-functions) which do not have this property, and then, the symbol hat will denote specifically 
the magic sheet. 

^This is a particular case of a general convention to avoid the default cuts of a function [34]. See [34] for more details 
about the shift operations. 

®In this paper, sometimes what we call the "analyticity" strip is in fact the "meromorphicity" strip since the corre- 
sponding function might have poles in that strip. 
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(a) The functions x(u) and l/x{u). (b) The functions x{u) and l/x{u). 



Figure 2: Riemann sheets of the multivalued function x{u). The subfigure 2a represents the 
function u i— > x{u) (upper sheet) and u i-> l/x{u) (lower sheet), whereas the subfigure 2b represents 
the function u f-)- x(u) (upper sheet) and u i— > l/x{u) (lower sheet). When u is real, x^^^i is defined by 
an analytic continuation from the real axis, along the red path which avoids the cut Zq, whereas i^^^l 
is defined by a continuation along the green path which avoids the cut Zq. We see that a;[+^l = £^+^1 
whereas a:'"^] = l/a;[~^l. 



class of functions analytic in the interval ™ < Im u < " is denoted An ■ Usually it means that 
we have a pair of Z-cuts or Z-cuts at the boundaries of the corresponding analyticity strip. 



In addition to the shift operation we will use the fused product for a function: 



(2.10) 



This definition comes from the g-number notation: [s]d = , where D — e^^". So, e.g. [2]d/ 

(D + D-i)/ = /+ + /- whereas /I^Id = /+/-. 

Finally, for the functions with cuts we define the discontinuity on the cut as follows 



disc (/) = f{u + iO) - f{u - iQ) . 
In this paper we use the following Baxter- type polynomials and "Z-polynomials"^ 



M 



Uj) , B 



(±) 





For convenience we chose a normalization which assures the following simple relations: 

^(±)^(±) ^ (_l)Mg± _ 



(2.11) 



(2.12) 



(2.13) 



The roots Uj of the Baxter polynomials in (2.12) are called "Bethe roots". They are related to Xj 
by Xj = x{uj). Their number and positions completely define different states of the s[(2) sector of the 
theory to which we restrict our analysis in this paper. 



'^i.e., having a fixed Z-cut on their Riemann surface, apart from zeros 
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In the integral equations that we wiU write we use the following kernels: 



JC{u) = -^ ' '*C-H = '*C("+f)-'*C(^-f)' (2-14) 



Z{u,v)^~^y^^^J^— , Zsiu,v)^Z{u,v-f)-Ziu,v+f), (2.15) 
VKzM^f + V(^I^"1-;^-^) . (2.16) 



2tt 27r ^ V 27r(n + 1) / ' 

n=0 ^ \ ' 1/ 

The convolution is denoted, for arbitrary kernels K and function /, by 

/oo 1*00 
K{u,v)f{v)dv, {^*f){u) = 4 K{u,v)f{v)dv, 
-00 'J —00 

K{u,v)f{v)dv, {K*f){u)= / +/ K{u,v)fiv)dv. (2.17) 

3 Fundamental properties of T-functions 

The purpose of this section is to find T-functions in a gauge with good analyticity properties, different 
for each of three bands of the T-hook (see figure 1), and then to relate the three bands by transition 
functions (gauge transformations). We describe and motivate a new structure in AdS/CFT Y- and T- 
systems which we call Z4 symmetry. It is tightly related to the analytic properties of T-functions which 
we consider first. We will also introduce a concept of the "magic" sheet. Although some analyticity 
requirements of this section look a bit speculative, they are verified in appendix C from the TBA form 
of the Y-systeni. 

It is important to stress that one cannot find a global gauge with good analyticity properties for 
all T-functions in the whole T-hook. Therefor we rather describe three different, "physical" gauges 
with a good analyticity for, respectively, the right, left and upper bands of the T-hook, which are then 
related by the transitional gauge functions. Then we show that the analyticity requirements to the 
functions of these gauge transformations appear to be strong enough^ to fix the system of FiNLIE and 
the corresponding physical solutions of the Y-system. 

3.1 Analyticity strips 

Whereas the Y-system itself follows to a large extent from the symmetry of the integrable theory 
the analytic properties of the Y-functions only partially are constrained by the consistency with the 
Y-system. However, they cannot be easily obtained from general properties of the theory. For the 
integrable theories for which the integrable lattice formulation is known the analytic properties can be 
read off from the underlying microscopic Bethe ansatz equations. However, in the present case such a 
formulation is unknown. Nevertheless, one can make reasonable assumptions based on the asymptotic 
solution of the Y-system proposed for an arbitrary state in [7] and from the integral TBA form of the 
Y-system for a subclass of sl(2) operators written in [15, 22]. 

The analyticity properties of Y-functions which are expected to be satisfied by the physical solu- 
tions of the Y-system are written in table 1, in the notations of the previous section. These Y-functions 



^together with some simple assumptions about the behavior at m — > 00 and about the structure of poles and zeros 
for a given state. 
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Ya,±l e Aa-1 , a > 1 

Yi^±s e A\s\-i , \s\>l 

Ya.s — Ya.s (Y-functions are real) 

Table 1: Analyticity properties of the Y-functions of the Y-system for AdS/CFT 

are considered on the mirror sheet which means that we have chosen the long Z-cuts [16]. Moreover, 
Yi ±1, 12, ±2 have a Z-cnt on the real axis and are related on this cut by 

i^2,±2(w + i0) = l/Yi,±i(M-iO) , ueZo. (3.1) 

The last property was discovered from the TBA form of the Y-system in [15] but it will be considered 
here as one of the basic analyticity assumptions for our construction of FiNLIE. It means that i^2,±2 
and Yi^±i are not independent functions but rather the analytic continuation of one into another 
through the ZQ-cut. 

The above list of analytic properties for Y-functions also imposes certain analytic properties for T- 
functions. Most of them can be inferred from the definition of T's (2.2), the analyticity of Y-functions, 
and the consequences of Hirota equation (2.5). Recall that the T-functions are not uniquely defined 
given the Y-functions, due to the gauge freedom (2.3). This makes the statement about the analyticity 
of T's a bit more complicated and gauge dependent. A natural assumption, which will be later 
demonstrated, is that, in a certain gauge which we denote by the bold font, Ta^s has good analyticity 
properties for the upper band, whether as in another gauge, denoted by the "blackboard" font, To^^ 
has good analyticity properties for the right band of the T-hook.^ These analyticity properties are 
explicitly written in table 2. 

Let us notice that there exists still some residual gauge freedom which will not spoil the analyticity 
properties. We will fix the remaining gauge ambiguity in the definition of ^ and ^ later on. In 
particular, we choose 

To,. = 1. (3.2) 

Now, with this knowledge of the analyticity strips, we will proceed with fixing the details of the 
analytic structure of T-functions in each strip. 

3.2 Group theoretical constraints 

Let us remind that there is some evidence [17, 30] for interpretation of the T-functions as transfer 
matrices. In particular, at strong coupling the T-functions are reduced to the characters of classical 
monodromy matrix, in such a way that each pair of indices (a, s) corresponds to a unitary representa- 
tion of the symmetry group PSU(2, 2|4) having a rectangular ax s Young tableau. This interpretation 
itself singles out a particular gauge. Indeed, according to [17] the pairs {n, 2) and (2, n), or (n, —2) and 
(2, — n), correspond for any n > 2 to the same typical representation and should have equal characters. 
For the transfer matrices this also suggests a natural gauge where 

T„,2 = T2,„ , T„ _2 = T2,-„ , 71 > 2 . (3.3) 
This is indeed the gauge in which the Wronskian solution for the T-hook was written in [18]. 

^For the left band the gauge is also diflferent from the upper and right bands but we will mostly consider the right 
band since the construction for the left band is similar. 
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if a > 


if ,s > a, 


Tq.O G Aa+l 




Ta,±l e Aa 




Ta,±2 e Aa-1 


Ta.s, Ta,s are real functions. 



Table 2: Analyticity properties of the T-function of the T-hook for AdS/CFT. 



Another constraint which we wih impose is Qg 1. In the rational (super)-spin chains Qg is the 
Baxter Q- function on the final step of the Backlund procedure [35, 36], where it corresponds to the 
trivial gl(0|0) sub-algebra and therefore it is naturally represented by a constant that does not depend 
on the spectral parameter. The same observation holds when Q-functions are constructed from the 
first principles, i.e. as eigenvalues of Q-operators [37-39]. Even though the construction of Q- and 
T- operators in AdS/CFT is still in a preliminary stage (see [30] for the first nontrivial quantum 
computation with T-operators) , we believe that will play the same role, i.e. it will be an identity 
operator. 

The quantum super-determinant should be also set to 1. In [18] the quantum determinant was 
identified with the following combination of the Q-functions: J^-Jl ■ Setting it to 1 and recalling that 

To,s = '^^ gives the following "unimodularity" constraints^ 

-•-0,0 ^ -"-0,0' -"-o.s — -"-o^o ^ -"-o.o • \'^-^) 

This tells that Tq.q should be an i— periodic function^^. Due to the importance of this quantity, in 
what follows we introduce a special notation 

T=^/T^. (3.5) 

We expect T to be a kind of a physical gauge. In this gauge the T-functions exhibit all natural 
symmetries of the characters, and for the left-right wing symmetric states (in particular, the s[(2) 
or su(2) states; we will call them LR-symmetric from now on) it can be shown to be essentially 
unique when supplemented by a set of analyticity conditions summarized in (3.44) (see Appendix E.l). 
Moreover, we will show that in addition to the properties listed above there is one more very important 
symmetry of the T-gauge which we call Z4 symmetry. 

Finally, let us define T in terms of T. A natural candidate for that is the following, very particular 
gauge transformation (2.3) from one to another: 

Ta^s = i-ir^'+'^TaA^^'^+'^r'' . (3.6) 

This transformation does not change indeed the gauge invariant Y-functions. The relation (3.6) can be 
considered as a minimal gauge transformation from T (where Tg^s, for instance, has an infinite number 

^"see the eq.(4.6) and the beginning of sec. 5. 3 there 

^^A closer look reveals an intimate connection of To,o with the dressing kernel in TBA. This is similar to the case 
of relativistic sigma models with rational S-matrices, where the scalar factor of R-matrix needed to set the quantum 
determinant to one is equal to the dressing factor of the factorized S-matrix in the integrable field theory with the same 
symmetry group [40]. 

^^We recall that since the original Y-system is defined on the mirror sheet, the corresponding Ta,s are defined as 
functions with long cuts. The periodicity condition is sensible to the choice of cuts. Periodicity of To,o can be also 
obtained from the TBA equations, see appendix C.2. 
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of branch points) to a gauge with the analyticity bands of table 2. Indeed, we see that To.s = 1- The 
property T2,±s € As-i is also obviously true since T2.a = T2.a- Then Ti.s = ±Ti.s/J^[-'^+''l seems to be 
a natural choice since both Ti ^ and Ti i should be regular functions^'^ with the same analyticity strip, 
according to table 2, and the factor will not spoil it. Moreover, as it is shown in appendix C.2, 
the property Ti^^ G As and (3.6) itself, for which we gave a motivation above, rigorously follow from 
the TBA equations. The overall sign (— l)°(''+i) does not affect the analyticity strips and can be 
chosen arbitrarily, by adjusting definition of the T-gauge. Our choice will eventually provide a simpler 
behavior of Q-functions at u ^ oo. 

3.3 Magic sheet 

Let us consider T-functions in the right band of the T-hook, i.e. for Ta^s with s > a. In this case the 
T-functions satisfy precisely the same Hirota equations as for the su(2) principal chiral model in [4]. 
Moreover, at least for the s[(2) states the Y-functions obey a similar asymptotics Yi^s — > s'^^l, u oo. 
This property is well seen in the ABA limit, see app. B.l, and also holds at finite size, see app. B.2. 
Therefore, to construct the T-functions with a good analyticity in this band we choose a gauge^^ T 
with 7i s — s, w — oo (i.e. it tends to the dimension of representation in this limit), preserving the 
analyticity strip As and the property 7o,s — 1- 

Following the logic of [4, 6] we find the general solution of Hirota equation satisfying all these 
properties in the right band of the T-hook in terms of the Q-functions^^ QitQitQ2,Q2 in ^ gauge 
where 

Q2-(3i = l, Qi--m + G(u), Q-2 = -Wi = -iu-G{u), (3.7) 

namely, 

= Q^''+0T'"'' , (3.8) 
where we introduced two mutually complex conjugate resolvents: 

G(u) = r dv^, lm(u) > 0, 

G{u) = +^ r dv(^, Im(«) < 0. ^ ' 

The function G (function G) is analytic everywhere above(below) the real axis and should be defined 
by analytic continuation in the lower(upper) half-plane. In general, one would expect G to have 
infinitely many Z-cuts Zq, Z-2, -2^-4, ■ • ■• Also notice that G^+'^l + G^^^l = p. In these notations we 
get 

To,! = 1 

Ti,. = .s + G[+^1 +G[-^1, (3.10) 
T2,s = (1 + G[-'+i1 - G["-i1)(1 + G[-"-il - G[-^+i1). 

In this way we write all T-functions, and thus Y-functions of the right band in terms of one single real 
function p. This result was first time presented at the talk [41]. Later an alternative way to reduce 
the right band to a finite set of functions was proposed in [42]. 

^^we frequently use the words "regular" and "analytic" in the paper for the functions without cuts but potentially 
with poles 

^^which will be later related to T-gauge by a gauge transformation conserving its analyticity properties 
^^in the notations of the paper [18] 
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Figure 3: Comparison of the mirror (on the left) and the magic sheets (on the right). 

71, s coincide on the real axis of both sheets. 



The function p can be fixed by gluing this right band of the T-hook to the rest of it, into the single 
T-system. A gauge-invariant way of doing this is to consider the following combination of Yi_i and 

^2 2 ■ 

A nice feature of this combination is that it can be expressed entirely in terms of T-functions from 
the right band: 



r= ; 1 ■ (3.12) 



Using (3.10) we find 

(l + Gl+21 -G)(H-G[-2] 



(3.13) 



(l + Gl+21 +G)(1 + G[-21 +G)' 
Due to the property (3.1) we should get 

r{u + iO) = \/r{u-iQ), u e Zq. (3.14) 

Since on the real axis the r.h.s. of (3.13) has branch points due to G and G, whereas G^+^l and 
are regular there, the relation (3.14) can be satisfied by imposing the following simple condition: 

Gl+o] = -G[^°l , ueZo. (3.15) 

This simple observation has some dramatic consequences. First, due to p = G[+°1 + G^-°^ we see that 
the density p has a finite support on the cut [—2(7, 217]. This allows us to clarify the definition of the 
resolvent as follows: 

2Tn J_2g u-v y~G{u) , lmii<0 

We clearly see the advantage of the choice of the Riemann sheet with only "short" Z-cuts Z„ over 
the "long" cuts Z„: the same function G{u), when being defined as a function with the short cuts, 
appears to be G(u) — a function with only one single cut! 

We can go further and convert any function with long Z^-cuts into a function with short cuts by 
the appropriate choice of the section of the Riemann surface (a priori of an infinite genus) . As we shall 
see, this strategy leads indeed to an essential simplification of analytic properties of various functions 
involved in the TBA equations. By this reason we call this section of the Riemann surface the "magic" 
sheet. We systematically denote various functions defined on that sheet by a hat over them leaving 
the same quantities on the mirror sheet without the hat. 



- 11 - 



3.4 Z4 symmetry 

In this subsection we study the properties of the magic version of T-functions in the right band of the 
T-hook which we denote by Ta,s- To produce this function we should simply systematically replace 
all G and G in (3.10) by G: 

to,. = 1 , 

Ti,. - s + G[+"1 -G[-"1, (3.17) 
h,s - /i.i ■ 

Now we present some important observations. Firstly, 7i,s coincide with 7i,s everywhere inside 
the strip As- This in particular means that 

ri+ ti- = r+ Ti- , s > 1, (3.18) 

thus the Hirota equation for Ti,. also implies the Hirota equation for 7i.s>i. However, for s = 1 this 
Hirota equation should be different since T^^t 1 i- s = 1 is precisely the point where the 

T-functions from the upper band start entering the Hirota equation. 

Secondly, although the initial parameterization (3.10) of Ta.s (and thus of Ta.s) makes sense^^ only 
for s > a, we can formally define Ta,s for any —00 < s < 00 by (3.17). Interestingly, for s = 1 we get 
then 

titiTi7i = ro,ir2,i + o, (3.19) 

which is consistent with 

Ti.o^O. (3.20) 

In fact, it is easy to see that Ta^s defined by (3.17) satisfy Hirota equation for any s if we define it for 
—00 < s < 1 as an analytic continuation in s.^^ 

Finally, the magic T-functions for the right band defined in this way satisfy 

fa.-s = {-irfa,s ■ (3.21) 

We call this new symmetry of the AdS/CFT T-functions the Z4-symmetry. In the appendix A we 
motivate its relation to the Z4 symmetry of the classical monodromy matrix of the string sigma model 
(inherited from Z4 symmetry in the super coset action construction). 

It is evident that T is related to T by a gauge transformation. Since in both gauges To.s = 7o,s = 1 

and T-functions are real we have to take in (2.3) ^4 = 1/gi, .93 = l/.g2, ^ — f = ^- Then, since 
in both gauges the analyticity strips are the same and given by table 2, h should be analytic in the 
upper half-plane. In addition, in appendix C.4 we show that h should be real on the magic sheet^^. 
Therefore T's can be found from the relation (3.17) as 

To,s = fo,s = 1 , ti,s = U+'^U-'^fi^s , t2,. = U+'+^^U+'-^^U-'+^^U-'-^'^f2,s, (3.22) 
from were it follows that Ta^s also obeys Z4 symmetry property given by (3.21). 



^^We remind here that the determinant expression (3.8) only holds when s > a. As explained in [18], the T-functions 
on the T-Hook are given by different determinants in each of the domains s > a, s < —a and a > |s|. 

^'^Let us point out that this analytic continuation of T-functions to negative s has nothing to do with the physical 
solution of Hirota equation on the T-hook at negative s. 

^*The function h on the magic sheet, denoted by h, is defined by equality h = h in the upper half-plane and taking 
all the Z-cuts being short. 
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Note that reality of h together with its analyticity in the upper half-plane means that h has only 
one single cut Zq, similarly to G. Therefore T-gauge inherits another "magic" property of the T-gauge 
— finite number of short Z-cuts. In particular, Ti has only two of them: Zs and Z^s- 

3.5 Z4 invariance of the upper band 

In this section we argue that the Z4 symmetry also leads to an additional symmetry transformation 
in the T-gauge having good analyticity properties in the upper band. In order to see it we have to 
pass to the magic sheet and then analytically continue the solution of Hirota equation in the variable 
a from a > \s\ to a < \s\, defining in this way Ta,s for an arbitrary integer a, positive or negative, 
providing that —2 < s < 2. Our claim, stemming from the Z4 symmetry, is that Ta.s defined in this 
way obeys the property similar to (3.21) of the right band: 

Ta,s = (-1)'T_,,, . (3.23) 

Let us demonstrate it. Consider a combination of Y-functions Yi 11^2.2 which reads in terms of T- 
functions as follows: 

Here we use that, by definition of the bold gauge, T2.3 = T3^2 and To,i = T^qq = T^qq- Next, from 
(3.1) we notice that ^1,112,2 gets inverted when passing through the cut Zq: 

Yi,i{u + iO)Y2,2{u + iO)^l/{Yi,i{u-iO)Y2,2{u-iO)) , \u\ > 2g . (3.25) 

If we simply substitute expression (3.24) into (3.25) we get^^ 



rp rp[+l-0] 

J-l.O -"-0,0 



-•-0,0 



(3.26) 



This relation takes a nice form in terms of T's. Employing the definition of the hatted functions we 
have Tq q^^*'' = Tq q^' (for |u| > 2g) from where we obtain on the magic sheet simply 

'i'(toTo,o ~ '^'1,0 ■ (3.27) 
On the other hand, Hirota equation on the magic sheet reads as follows: 

T(to'i'o,o — Ti,oT_i^o + ToaTo,_i , (3.28) 
which coincides with (3.27) if (3.23) is satisfied! Indeed, from (3.23) we have 

to,±i = (3.29) 

and 

ti,o = t_i^o. (3.30) 

To give some more evidence to these results let us consider a bit more complicated example. Take 
another expression built from Yi^i and I2.2 which is rewritten entirely in terms of the T-functions of 
the upper band of the T-hook (a > \s\) 

l + i^2,2 _ TJ2T2 2Ti,o 



1 + l/Fi,! T+iTr,iT3,2 

^®we use a natural notation j["±01 = f(^u -\- in/2 it iO) 



(3.31) 
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From the property (3.1) of F-functions it follows again that 



s{u + iO) ^l/s{u-iO), \u\>2g. 



(3.32) 



If we calculate s slightly bellow the mirror Zq cut we can represent it as a product of the following 
two ratios: 



and we see that (3.32) is satisfied assuming, as usual for the AdS/CFT integrability, that all branch 
points are of a square root type (in fact only Z-cuts are present). This is an additional argument for 
the hypothesis that To,i = . 

In appendix C the Z4-symmetry of the upper band is derived rigorously from the TBA equa- 
tions. However, from our point of view, the Z4-symmetry should be rather included into the list of 
fundamental properties of the AdS/CFT Y-system, out of which the TBA equations can be derived. 

3.6 Bethe equations 

In the study of integrable spin chains a very convenient way of writing the Bcthe equations is the 
Baxter equation supplemented with the condition that the eigenvalues of transfer matrices are analytic. 
Analyticity is an obvious fact from the very definition of the transfer matrix, but it is not immediately 
clear from the explicit expression for its eigenvalues in terms of the Bethe roots and is only true once 
the Bethe equations are satisfied. 

Provided that the interpretation of [20] is correct and the T-functions of AdS/CFT are indeed the 
transfer matrix eigenvalues one should expect that the auxiliary Bethe roots (carrying no momentum 
and energy) can be found by requiring some good analytic properties for the physical transfer matrices. 
In this paper we restrict ourselves to the sl(2) sector, which has no auxiliary roots, and thus the 
analyticity should be simpler. 

As it is shown in appendix E.l, one can make gauge transformations by means of i-periodic 
functions without Z-cuts preserving all the properties of T- and T-gauge mentioned before in this 
section. We claim that it is possible to show that the requirement of absence of poles in the T- 
functions in T- and T-gauges, at least in their analyticity strips, together with the requirement that 
To,o has the minimal possible number of zeroes in its analyticity strip, removes the residual gauge 
ambiguity. This fixes the T- and T- gauges completely up to an inessential constant factor. 

Now suppose that Tq.o has some number of zeroes in the strip —1/2 < Im (w) < 1/2. Since 
J- — y^^Tofi defines the gauge transformation between T- and T-gauges, the only possibility to avoid 
the appearance of branch points different of those of standard Z-cuts, is to have only the double zeroes 
in To,o, so that T has only simple zeroes. We denote these zeros as Uj and assume that there are M 
such zeroes. 

Comparison with the TBA equations shows that uj are nothing but the Bethe roots and thus we 
should also satisfy 




relations give 'it.2'^2.2 = T3,2Ti,2 and, if To,i = then TjiTr,i = Ti^aTi^o • Thus 

rr[-l + 0]rp[-l-0] 




(3.35) 
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Figure 4: Paths used for analytical continuation. To get the Bethe equations used in the TBA 
approach, we continue Yi o using the path 7 (on the left). Alternatively we can formulate Bethe 
equations using Ti 1 continued over and 7^ (on the right). 



where Y^q denotes the analytic continuation along the contour 7 defined on figure 4. 

Let us mention a curious observation. Since Ti_2 ~ Ti_2/^^, the absence of poles in Ti^2 is 
only possible if T1.2 has zeroes at positions Uj ± i/2. Assuming that T2,i does not have zeroes at 
Uj ± i/2, ¥2,2 = T2,i/Ti,2 should have poles at uj ± i/2. On the other hand, let us consider analytic 
continuation along the contour 7 of the Y-system equation at a = 1, s = 1. Using (3.1), one gets 



^/y2. 



On the l.h.s., the poles of the 2 zeroes of the Bethe roots are conveniently canceled with zeroes 
of (1 + Y^q), due to (3.35). Therefore we see that the exact Bethe equations (3.35) can be replaced by 
the condition that the r.h.s. should be regular at u — uj. We believe that it reduces to the condition 
of regularity of (Y^^)'' at it = uj. This is an interesting statement deserving a further study. 
In appendix E.2 we show that the Bethe equations (3.35) can be alternatively written as 



f7+ 
1,1 

^1,1 



(3.36) 



where the analytic continuation along the contours 7± is defined on the right of figure 4. 

Eq. (3.36) looks more natural than (3.35) since Ti^i has only one Z-cut in each of the half-planes 
of the magic sheet, therefore the analytic continuation along 7_|_ or 7_ is uniquely defined. On the 
contrary, Yi.o has infinitely many ^-cuts in the mirror sheet and one should impose an additional 
condition in (3.35) that the contour 7 goes only through the closest cut. In appendix E.2 we analyze 
(3.36) and represent it in terms of the quantities suitable for numerics. 

3.7 Expression for the energy 

The transfer matrices of integrable models form a commutative set of operators. Their expansion with 
respect to the spectral parameter produces a family of conserved charges, one of them being energy. 
Hence we expect the energy to be encoded into certain asymptotics of the T-functions. To see how it 
happens in AdS/CFT case let us first consider the TBA equation for Yi.1^2,2 given in [15] (the sum 
of eqs. (40-41) therein) 



M 



Yi,iY2^2 = n 



exp 



; Za * l0g(l + Yafi) 



a=l 



Xk = x{uk). 



(3.37) 
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The kernel Za is defined in section 2.2 and for u — )• oo + zO can be written as^ 



The product in the pre-exponent has a similar expansion: 



exp 



.fc=l 



Therefore the exact finite volume expression for I1.1I2.2 has the following large u expansion: 

iE 



i^i,i^2,2 = exp 



(3.38) 



(3.39) 



(3.40) 



where we used (2.9). Note that for all physical states P — Q 

Ti,o 



Now let us recall that Fi. 1^2,2 — ^tf^- Since Ti.o is regular on the real axis we get due to (3.1) 



log 



To,o(u 



iO) 



To^o(u- 



-21ogri,ir2,2 , \u\>2g. 



(3.41) 



This important property essentially defines the function Tq.o in terms of the product of Y-functions 
(see appendix E.3), and it will be often used in this paper. Since Tq^q is an z-periodic function, one 
gets: 



, To,o(u-|+iO) 

log 

^To,o(w+|-iO) 



-2- 



iE 



(3.42) 



(3.43) 



so that both Tq.o's are inside the analyticity strip. We can see that this expression implies 

E ^ \ hm udu log Tq.o ■ 
Thus To,o renders indeed the value of energy of the state when being expanded in u. 
3.8 Summary of properties of Ta,s and ^ 

In this subsection, for convenience of the reader we summarize the properties of the T-functions 
encountered above. 

First, the T-functions satisfy the following analyticity properties: 



^"the next terms contains the next "local" charges. For example the 1/n^ term is 81193 ,a(")+a ^ Note that also a 



"non-local" part ~ a is present. Here 5„ 



^ are "local" charges [32]. 
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Ta,0 G Aa+1 




T's are real functions analytic in the upper band: 


Ta,±l G A-a 
Ta,±2 e Aa-1 


(3.44a) 




T„,2 = T2,„ , n > 2 




T's satisfy "group-theoretical" properties: 


T„ -2 = T2 _„ , n>2 

-■-0,5 — -1-0,0 


(3.44b) 


Ta,s obey Z4 symmetry on the magic sheet: 




(3.44c) 


T's have no poles in the analyticity strip 




(3.44d) 


Tq^o has a minimal possible amount of zeroes 




(3.44e) 


The double zeroes of To,o are the Bethe roots Uj 




(3.44f) 


Then, the T-functions are related to the T-functions by a gauj 


le transformation involving 




Ta,s = (-1)"(-'+1)T„,,(J-['^+^1)'^-2^ 


J' = \/To,o ■ 


(3.45) 


They satisfy the properties: 






T's are real functions analytic in the 
right (left) band: 


'^a,±s = 1 

T2.±s e As-i 


(3.46a) 


Ti^s has only two magic cuts Zg and Z-s 




(3.46b) 


Ta.s obey Z4 symmetry on the magic sheet: 


Ta.s = 


(3.46c) 


T's have no poles in the analyticity strip 




(3.46d) 



In this paper, restricted to the states of s/2 sector, the symmetry between the right and left wing 
implies an extra relation T!a,s = Ta,-s- 

All the listed properties of T-functions can be derived from the TBA equations, as we show 
in appendix C. However, in this paper we take a different point of view and consider them to be 
fundamental. Below we fix the solution using only these properties and taking some extra information 
about the large u behavior of various functions from the large volume asymptotic solution. 

In appendix E we show that the conditions (3.44) assure that the T-gauge is unique up to a 
normalization constant. 

4 Wronskian solution 

In section 3.3 we managed to express all T-functions Ta,s in the right band of the T-hook in terms of 
a single function — the spectral density p. In this section we will find a similar representation for the 
upper band of the T-hook by using the so-called Wronskian determinant solution of Hirota equation 
in the infinite band —2 < s < 2, —00 < a < 00. The right (left) and upper bands of the T-hook will 
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Figure 5: Infinite band i?'"-* for tlie Hirota equation. 



be represented by 2 x 2 and 4x4 Wronskian determinants, respectively. In the next section, the full 
finite set of equations, FiNLIE, will be found by gluing these three bands together into the full T-hook. 
Each of these steps has to be done by respecting the structure of analyticity strips of T-functions. 

The Wronskian solution allows to parameterize the infinite set of T-functions satisfying Hirota 
equation in a band, in terms of a finite number of Q-functions. It is thus important to understand 
the analyticity properties of the underlying Q-functions in virtue of the analyticity conditions and the 
symmetries, such as Z4 symmetry, established in the previous section. The representation in terms 
of Q-functions will be the base for construction of the FiNLIE system for AdS5/CFT4 formulated in 
sections 5,6 and solved then numerically in section 7. 

4.1 General Wronskian solution 

Here we describe the Wronskian solution for Hirota equation (1.1) on an arbitrary infinite band of 
width n shown in figure 5. Let us denote this band as B^"\ This kind of Wronskian solutions was 
used in [35] for the analysis of the su(n) quantum spin chains and in [6] for the solution of the su(n) 
principal chiral field (PCF) model in a finite volume. At this point, we do not specify any particular 
boundary conditions at the end (s = 0) of the band (we need three semi-infinite bands, one 5(4) and 
two B'-^\ for the construction of the T-hook). 

A possible basis of the Wronskian ansatz is a set 2n -f 2 Q-functions: 90, qi and p0, pi where 
i = 1, 2, . . . , n. To get rid of numerous indices it is very convenient to use the formalism of exterior 
forms. Namely, introducing an auxiliary basis of n vectors e* we define the 1-forms: 

q= qte^ , p= pie\ 

Then, introducing the k-forms^^ 

q[k-i] /\ q[k-d.] _ _ , 

= [k-2] [fc-4] ' ' ' [2-fc] ' 5(0) = 90 ' (4-1) 

% % ■■■% 

(and similar definitions for p(fc)), we write the general Wronskian solution for the _B(") band [35] simply 
as follows 

Ta,s = A pi;!) , (4.2) 

where we identify the exterior n-form with a scalar, A • • • A e" = 1. One can easily check that for 
arbitrary q's and p's, T^^s defined in this way satisfy indeed Hirota equation with Ta.s = outside 
the band i?'"'. Notice however, that there is a certain freedom in choosing different sets of gg, qi and 

^^In spite of this similarity in the Wronskian representations, the analytic properties of Q-functions are totally different 
in AdS/CFT and PCF models, especially due to the difference in position of the momentum carrying nodes in the band. 
1(.k) = h.Hi-ik^'^ A...Ae*fe. 
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-■ ,^-6 — 6 — 6— 



Figure 6: The structure of Z4 symmetric band B^-^^ obtained by analytic continuation from the 
right band of the Y-system. T- and Y-functions to the right of the dashed diagonal and T-functions 
at the dashed diagonal are the ones of the AdS/CFT T-hook and are given by the Wronskian 2x2 
determinants. T-functions to the left of the diagonal are the analytic continuation of these Wronskians. 

P0, Pi for the same solution of Hirota equation. Indeed, the construction possesses the following sl{n) 
symmetry: 



{detH = 1) 



(4.3) 



which leaves T^, invariant for an arbitrary non-degenerate i-periodic matrix H{u). There are also 
two scalar symmetries: 



P'i. 
Pi 



Fp^, 



Pit, C^'pdi 
Pdi Fpq, 



F-'q^ 



(4.4a) 
(4.4b) 



with arbitrary ^-periodic functions C and F. 

Eqs.(4.3) and (4.4) form the complete set of transformations of the Wronskian solution which leave 
T's invariant. Indeed, qi,q2, . . . q„ should be n independent solutions of the Baxter equation (see for 
instance [35]). Any other solution of the Baxter equation should be a linear combination of these n 
solutions with coefficients being i-periodic functions, i.e. it should be related to those n solutions by 
a combination of transformations (4.3) and (4.4). The same is true for pi,p2 ■ ■ ■ ,Pn- We have also an 
additional freedom in rescaling p0 and (70 which explains why there are two scalar symmetries (4.4) 
and not one. 



4.2 Right band 

In this subsection, to demonstrate to the reader the method of Wronskians, we rewrite the results of 
section 3 for the right band in Wronskian notations. Also we will discuss rigorously the analyticity 
properties of Q-functions of the right band previously assumed in section 3 as being just natural. All 
this will be very helpful when we generalize these methods, in the next couple of subsections and in 
appendix D, to the considerably more complicated case of the upper band. 

Using Eq. (4.1) we can write the T-functions in the right band in terms of 6, so far arbitrary 
functions (70, 



91, 92, Pi, P2^^ 



Il2,s-P0 qi 



i(2) 



where 



9(2) 



q+ hq 



P(2) 



p^ Ap 



(4.5) 



(4.6) 



90 P0 

This is the most general solution of Hirota equation in the infinite horizontal B^'^^ in (a, s) plane, 
without any analyticity assumptions. Let us show how Z4 symmetry constrains the form of this 
solution and its analyticity properties. 



^In the notations of eq.(4.11) of [18] we liave • 



■ Q0, "ji = Qi, 92 = Q2, P0 = Qt2, pi = 02/-^, P2 = Qi/J^- 
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• From Ti.o — 0, which is a consequence of Z4 property (3.46c), we conclude that the forms fi and 
q are hnearly dependent: p — aq, where a is a function of the spectral parameter. 

• Ti l — ~Ti^_i implies that a+ ~ , which means that a is i-periodic. It is therefore possible 
to absorb a into p and q using (4.4b) with F = ^fa. Hence we put a = 1 in what follows. 

• From Tq.s = 1 we conclude that q^ = q^ and — ggTij. 

From all these properties we get the solution which depends only on qi and q2'- 

to,s = 1 , Ti.s = g'+^'i A , f 2,. - il+^lfl^l . (4.7) 

The solution (4.7) literally coincides with the parameterization of 5u(2) XXX spin chain transfer 
matrices in terms of the Baxter Q-functions^^ [35] . This implies that the following system of Baxter 
equations should be satisfied: 

1,1 = q+T^+;-'^ - rt'+li , (4.8) 

which is easy to check explicitly by substitution of (4.7). 

The equations (4.8) are very useful for the analysis of analytic properties of Q-functions induced 
by the analyticity of T-functions. Using them we can prove that the analyticity of Ti^^. inside the strip 
As and of Ti_i everywhere except Z±i, implies the existence of such a symmetry transformation H 
from (4.3) that qi and 92 are analytic everywhere except a single cut on the real axis. In particular, 
this shows that all Ti^s have only two magic cuts. 

Here is the idea of the proof. From the analyticity properties of T-functions and (4.8) it follows: 

f ^^=^1 disc g[2'-4] ^ Ylr-i] disc <j[-2l _ 'f^--^\ disc gl-^l ^ r>2. (4.9) 

In other words, the discontinuities of g's on the cuts Z_2 and Z_4 define all discontinuities on 
r > 0. 

We see that if a solution is regular at Z_2 and 2^-4, it is then automatically regular in the upper 
half-plane. The conjugate of equation (4.9) expressing the discontinuities of g's on Z_2r in terms of 
the discontinuities on Z2, Z^ implies the analyticity in the lower half-plane. It is enough to find such 
a symmetry transformation (4.3) that 

disc (H/ )qJ.^^"l + H,^ disc (q]"^"l) = 0, n 1, 2 , (4.10) 

which represents a Riemann-Hilbcrt problem. Assuming that it has a solution we can prove the above- 
mentioned statement. In order to argue that the solution exists let us recast Eq. (4.10) into a linear 
integral equation 

Hi = VI + ^ coth(^(^. - v))i [iJf disc {Al){A-')i] , A- = (^-^^^ , (4.11) 

where Vf are some z-periodic functions without cuts. It is convenient first to relax the condition 
detH — 1. Then (4.11) is a linear equation for four independent matrix elements of H. Suppose we 
found a solution of (4.11) which means that there is indeed a linear combinations qi of the initial qi 

^^up to relabeling Ta.s — ^ Ta, a-i- Of course, in our case Q-functions are not polynomials. 
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Figure 7: The structure of the Z4 symmetric band B^^^ obtained by analytic continuation 
from the upper band of the Y-system. 

with periodic coefficients without cuts Z_2 and Z_4. As it is easy to see, also satisfy (4.8) and, 
as a consequence, (4.9) which then implies that qi could have only one single cut Zq. We should 
remember, however, that the expression for Ti ^ in terms of qi will also contain detH in denominator. 
This denominator is easy to get rid of since it is a periodic function and cannot contain any cuts. 
Indeed, 

and so we can absorb detH into keeping the nice analyticity properties of gl+^l . 

It is important to mention that the explicitly found in [18] large volume asymptotics of Q-functions 
has indeed only one cut. 

We still have some freedom in the _ff-transformations for Hi^ being regular i-periodic functions. 
A part of this freedom can be used to make qi real and (j2 purely imaginary. Both functions have only 
one cut and thus can be very efficiently written in terms of a spectral representation (i.e. from the 
discontinuities on the cuts). 

We finish the discussion of the right band by relating the q's to the quantities introduced in 
section 3. One can easily identify qi = (— iu + G)h and q2 = h, where G was defined in (3.16). 

4.3 Upper band 

For the upper band we can perform a similar procedure, but now in the T-gauge. The following 10 
functions can be used as a basis: qg, pg and qi,Pi, i — 1, . . . , 4.^^ 

The multi-indexed Q-functions qij, q^fe and ciijki can be determined from the Pliicker determinant 
formulae (4.1) in the vertical i?'-^' band of the T-hook:^^ 

q+ A q^ q++ A q A q qI+^1 A q+ A q^ A ql^^l / . , „x 

q(2) = — - — , q(3) = , q(4) = ++ — — ■ (4.13) 

q0 %% q0 q0q0 



^^In the notations of [18] they arc q0 = Q12, P0 = Q34, qi = Qi2i' Pi = Q34i- 

^^There exists a Wronskian parameterization of the general solution in the full T-hook [18, 28] and it can be also 
represented in terms of exterior forms (we postpone it to future publications) . But in this paper we prefer to build our 
physical solution from three Z4-symmetric bands. 
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Wronskian solution ol the vertical band is then written in an extremely compact form^^: 

Ta,s=4l[^Ap[Z[)- (4-14) 

The Wronskian ansatz (4.14) gives the formal general solution of Hirota equation in the vertical 
B^^^ shown in figure 7. Further restrictions on the underlying q's are needed to ensure that the Z4 
properties are satisfied, as well as the other analyticity properties listed in section 3.8. Moreover, for 
simplicity we consider the LR-symmetric states only^* and hence an additional, LR wing exchange 
symmetry is imposed; 

t„,_, = t,,, . (4.15) 



In appendix D we study all these conditions in detail. For the moment let us notice that from 

n[+a]pTp[-a; 
1,1 1,1 



(3.45), (4. 7) we have fa,2 = T2,a = = t^^^T^ °1 and ta,2 = ta,-2- Therefore we can set 



qfl = P(4) = P0 = q(4) = TTi,! , (4.16) 

by making appropriate transformations (4.4). 

After that we still have a residual s((4) symmetry (4.3). As we prove in appendix D.8, this 
symmetry can be partially used to choose q-s and p-s satisfying the following relations: 

qi23 = qi = ^Pl34 = — P3, qi24 = q2 = P234 = P4, 

qi34 = q3 = -Pl23 = -Pl, q234 = q4 = Pl24 = P2- (4.17) 

As a consequence of (4.17) we also get 

qi2 = P34, qi3 = Pi3, qi4 = -P23, q23 = -Pl4, q24 = P24, q34 = P12 ■ (4.18) 

With these identifications, we obtain from (4.14) an explicitly real and LR-symmetric parameterization 
[18]: 

ta.il = q^+'^lqi"'^' + q^^^'q^"' + qli+'^'q^'^' + qi""°'4""' , (4-19) 

-■-0,0 - qi2 qi2 +q34 q34 -qi4 qi4 -q23 q23 -qi3 q24 -q24 qi3 ■ 

One may wonder whether the Wronskian solution (4.14) possesses a finite cut structure for the 
Q-functions, as it was the case for the right band. We performed a detailed analysis of this question 
in appendix D and came to the conclusion that unfortunately at least some of Q-functions should 
have infinite number of cuts. However, we were able to show that there is a choice of Q-functions 
such that q(fe) is analytic in the upper half-plane above .^-i+|fc-2| ^-nd P(fc) is analytic in the lower 
half-plane below 2^i_|fe_2|i where the Z-cuts are absent. This analyticity condition fixes a part of the 
s[(4) symmetry. It can be shown that the transformations (4.3, 4.4a, 4.4b) used to enforce (4.17) do 
not spoil this analyticity condition, so that we can impose that (4.17) holds for Q-functions such that 
q(fe) is analytic in the upper half-plane above Z_i+|;j„2|- In the rest of the paper we stick to this 
analytic choice of the Q-functions. 

A nice property of the relations (4.17), especially important for the numerical applications, is 
that for the large volume L the terms in the Wronskian formulas (4.19) are well distinguished by their 
magnitude: the first two terms in Ta.±i are of the order one whether as the other two are exponentially 



^^In this section we systematically remove hats over q-s and p-s for convenience, however we consider them as functions 

with short cuts. 
28n 



'This includes in particular all sl{2) and su{2) states, and many more. 
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small w.r.t. the length L. Similarly, the first term in Ta_o is of the order 1 whether as the other are 
exponentially suppressed. 

Let us show that it is enough to know qi,q2,qi2 to restore unambiguously all the q-functions 
of the upper band. For that we use the Pliicker relations which follow from (4.14) and are explicitly 
given in [18, 20, 36]: 

q0q*j = q^qj - '^t^l^ - (4-20) 

q^jfcqi = q.+ qrfc - q^q^ • (4-21) 
Reconstruction of the rest of q functions goes as follows: 

• From (4.20) with ij = 12 one finds 

q« = ^^^^ ^^^^ . (4.22) 

qi2 

• Relations (4.21) for ijk = 123, 124, 213 and 214 can be explicitly written in view of (4.17) as 

q? = qjjqis - qKjqr2: (4-23) 

qiq2 = flt2^u - <^U<^12 = <^t2^2Z - q23qr2' (4-24) 

q2 = q^2q24-q24qr2- (4-25) 

We can unambiguously and explicitly define qi3, qi4, q23 and q24 through qi, q2, qi2 using their 
regularity above by inverting the linear difference operators in these Pliicker relations. Note 
that qi4 = q23. 

• Similarly, the relations (4.20) for ij = 13 and 14 define qs and q4. 

• Finally, the relation (4.21) for ijk — 134 fixes q34. 

Let us summarize the analyticity properties of the Q-functions introduced in the previous section: 
qi,q2 are regular above Zq and qi2 is regular above Z_i. A simple inspection of relations (4.20) and 
(4.21) shows that this automatically ensures the correct analyticity of all the rest of q-functions, and, 
consequently, of all T-functions of the upper band: q0 is regular above Zi, q^ are regular above Zq 
and are regular above 



5 Finite set of equations 

In the previous section we managed to parameterize the T-functions of all three bands of the T- 
hook by Wronskian determinants of the Q-functions and establish their symmetry and analyticity 
properties. Here we will derive yet missing equations of the system FiNLIE gluing all three bands 
together into the single T-hook by means of a few transition functions and constraining these functions 
by their analyticity properties. The method will not appeal to the TEA equations but rather will be 
based on the properties listed in section 3.8, from which the TEA equations also follow (as is shown 
in appendix C). They look more fundamental and simple then the set of analyticity properties for 
the construction of TEA from the Y-system given in [11] but their equivalence is demonstrated in 
appendix C. 
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5.1 Density parameterization of the upper band 

In section 4.3 we explained how all the T-functions of the upper band can be written in terms of only 
three independent functions; qi,q2 and qi2. Let us now introduce a suitable parameterization for 
them: 

m^uf+r, ci2^uf+rw, qi2 = /'Q, (5.1) 

where Q — 11^=1 ~ ^i) simply a polynomial of degree M containing all zeros of qi2 and 
normalized so that Q = u*^ + . . ., whereas /, U, W are nontrivial functions. We see that this is indeed 
just a parameterization which does not change the number of independent functions, and the roots 
of the polynomial Q encode (but are not equal to) the Bethe roots of an s[(2) state. One can think 
about this parameterization as being a gauge transformation to a new gauge 3^a,s = 9(2^i) ^ 
defined by 

f + J- f++ j — 

q0 = U+U-f++f—qti, q = Uf+f-q, q(2) ^ .fq{2), 1(3) = -^j—q{3), = jj+jj- g(4) (5-2) 

so that: 

qi^l, q2^W, qi2 = Q, qi23 = ■ (5.3) 
The ^-functions and the T-functions are related as follows 

Ta,s = ^a,. 7'^"+'' (t/[+"li7[~°iy''" . (5.4) 

In particular, one can see that the LR symmetry of Ta,s implies that 

= (C/[+''lc7[~'^l)2=^a,i . (5.5) 
The Wronskian representation (4.19) becomes: 

X,i = VKl-^'l + W^I+^'l + + , (5.6) 

^a,0 - 912 1l2 +934 934 " 9l4 9l4 "923 ^23 ~ 9l3 ^24 ~ 924 9l3 • [^■') 

Knowing the analytic structure of the functions (5.1) we can parameterize them in terms of spectral 
densities with the support on M. Not only it is conceptually important for the reformulation of the 
whole Y-system as a finite Riemann-Hilbert problem but it is also very convenient for the numerics. 

First, as we see from (5.1) = ^ and thus it should be regular in the upper half-plane. This 
allows us to introduce the spectral representation for W. Defining a real function 

P2 = VK[+°1 + VK[-"1 , (5.8) 

for the states with two symmetric magnons we can write 

W = -iu + JC*p2 , lmw>0, (5.9) 
W^+iu-IC*p2 , lmw<0. (5.10) 

For the states with more than two magnons the linear polynomial —iu in W should be replaced by a 
polynomial of degree M — 1 (see appendix B for more details). 



In the large volume limit Q coincides with the Baxter polynomial Q{u) = njiiC" ~ ^j)- Recall that To,o = 
C|l2qi2 -I- subleading terms . The zeros Uj of qi2 are chosen in such a way that To,o has zeros at the positions of Bethe 
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This function W used for parameterization of the upper band is a direct analog of the function 
Qi parameterizing the right band and defined before in eqs.(3.7),(3.9) as 

Qi~—iu + JC*p. (5-11) 

In the large volume limit, the densities p and p2 become semi-circle distributions with a finite support 
on Zq. At a finite volume p still has the same finite support whereas p2 becomes non-zero everywhere 
on the real axis. However, if we define p2 = P2 +93^*''94 +94^'^' 93 °' then the Z4 constraint J^o.i — 
implies'^'^ that p2 has a finite support, in complete analogy with p. In terms of this density W is 
parameterized as follows: 

W = -iu + IC*p2-IC* (4^°' 9!"°' + , Imw > . (5.12) 

Contribution from (73 and (74 vanishes for large volume and thus the parameterization (5.12) is very 
handy for numerics. 

Similar spectral representation can be written for / and U and thus we conclude that the spectral 
problem reduces to a problem of finding a few densities as well as a few additional parameters, such 
as the Bethe roots. In the next subsection we show however that the function / can be explicitly 
excluded from the final FiNLIE. The further restrictions on the densities are due to the symmetries, 
analyticity properties and the condition that the Hirota equation is satisfied not only inside the bands 
but also in all nodes of the T-hook. In other words one should sew the three bands of the T-hook 
together to close FiNLIE. 

5.2 Closing the system of FiNLIE 

In this subsection we show how the system of equations can be closed by sewing the T-hook from three 
bands and imposing the analyticity properties discussed earlier. 

5.2.1 Equation for / 

The regularity of the Q-functions in a half-plane allowed us to write a spectral representation of W in 
terms of a single density p2 (5.12). Here we will exploit the analyticity properties of U and /, which 
follow from their relation to the Q-functions (5.1), to write the equations for fermionic Y-functions 
Yi ±1, l2,±2- Namely, we will use the fact that the functions U,f~,h are regular in the upper half- 
plane. Moreover, for a sufficiently large L (or sufficiently small coupling constant) these functions do 
not have poles and zeroes in the upper half-plane and behave as a power of u at m — > 00, as one can 
see from the asymptotic limit (see appendix B). 

Considering a function (^j^^ we see from (5.4) that it can be written through T-functions of the 
upper band and the fermionic Y-functions as follows: 

l/v 'Tuo.^o:o'yuiY2a^o:o' ^ 

Moreover, due to the analyticity properties of / this function is analytic in the upper half-plane and 
goes to 1 at infinity there. Also we assume (and our numerics seems to confirm it) that it has neither 
zeros nor poles there. Thus one can construct a spectral representation for log B in the upper half-plane 

\ogB = JC*pb , lmu>0, (5.14) 
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One can define p2{u) = lime_>+o ^e,iiu). For u S Zq the limit is zero since it is equal to ,^0.1 which must vanish 



due to Z4. From (5.6) we get p2 = P2 + 53^"' 94 + si^^'^a 
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from its real part on the real axis 



log + , \v\ < 2g 



Pb(v)=logB(v + iO)B(v-iO)=< --0^-0,0^1.1-2,2 ^ (5^5) 

log^^f^ , \v\>2g 

V. --o, 0-^0,0 

Alternatively, one can also reconstruct B 

logB = /C*% , lmu>0 (5.16) 

from its imaginary part 

= log ^^^^ ^ = < ^; / \2 . (5.17) 



B(.-^0) logfif^;,^^) , H>2g 

^_ --0,0 \^11 ^22 / 



Notice that in both representations (5. 14), (5. 16) the term logB(t; — iO) does not contribute in the 
r.h.s. We also used the fact that due to (3.1) log(yi^il2.2) is real between the branch points zt2g and 
is purely imaginary outside that interval. 

The former representation has an advantage with respect to the latter one since Pb{v) tends to 
zero faster than rji, when v — > ±00. This is because at large B is asymptotically a phase^^ as can be 
seen from (5.13). 

Since B is a simple combination of /'s, one can determine / from the following finite difference 
equation: 

log(r)2 ~ log(/+)2 = logB = /C * . (5.18) 
Keeping in mind that / should be regular in the upper half-plane we write the solution as follows 

00 

log/2^^/C[2"-il*Pfc. (5.19) 

71=1 

The sum in the r.h.s. is divergent, however the divergence is just an infinite constant which can be 
regularized, as in (2.16), so that 

log/2 = *+*Pb, (5.20) 

where the kernel ^ is defined in (2.16). Note that in principle one can add a constant to the r.h.s. 
of the previous equation. This however would change only the normalization of / and as a result the 
normalization of T. We fix the normalization by (5.20). 

5.2.2 Equations for l^i. 1,12,2 and U 

Equation for Yi^il2,2: To derive an equation for Yi,il2,2 we consider instead of (5.13) 

Again we will use the fact that due to (3.1) log(yi, 1^2,2) is real between the branch points ±2g and 
is purely imaginary outside that interval. Together with the square root the product Yi.iK2,2 drops 

^-"^In (5.13), the factor 5^i,o is real, Y11Y22 is a phase as soon as u > 2g and »9i,o/.3q~o approaehes 1 at large values of 

u. 



- 26 - 



out from the imaginary part Im (log B) and therefore does not appear in the r.h.s. of the spectral 
representation 

logB = /C * 77h, Im(M) > , (5.22) 

where 



^, = log — = <^ ~- -0.0 (5 23) 

Inserting the definition of B (5.13,5.21) into the l.h.s of (5.22) and shifting the contours of integration 
in the integrals involving log we obtain an equation for Y11Y22 expressing it through and 
5o,o and the polynomial Q{u) — YijLii'^ ~ Uj) encoding the Bethe roots of the state: 

log Yi, 1^2.2 = log + Z+a * log gTg^ ^ ^1 * "g^ ' '^^■2'^) 

where the kernels are defined in section 2.2. This equation can be also obtained from the TBA 
equations as is shown in appendix C.2. The derivation from TBA is an important check of our basic 
analyticity assumptions in section 3. 

Equation for 11,1/12.2: Similarly to li.il2,2, we can derive a separate equation for the ratio 
li, 1/12,2- For that we construct a combination of T-functions which contains only U and / and 
no conjugate quantities, namely 

Ti,o(Tr.i)^ _ / U ^1,0(^1:1)^ 

To-o(T2,i)2 =^o:o(=^2,i)^ ■ ^ ^ ^ 

Next, writing the ratio of the Y- functions in terms of T, 

' (Tl,2)^ 



Y2.2 



Ti,o(Tr,i)^ 
To:o(T2,i)2 



and using that, according to (3.45) and (3.22), ;^/Ti,2 — ^^7^,1/71,2 we construct from here the 
quantity 

l^c/[2i f^h ) " Y2,2 ^1,0 l^ri,2 ^i:J ^^^ ^^ 

which is written completely in terms of functions analytic in the upper half-plane and approaching 1 
at infinity. Its spectral representation is again straightforward: 

log C = /C * , Ini(w) > (5.28) 

with the spectral density given by 

Here we use the reality of T-functions as well as the reality of log Fi, 1/12,2 so that the Y-functions 
again drop out from the density. Furthermore, we can get rid of the shift in arguments of the functions 
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in the definition of the density by shifting the integration contours in (5.28). As before, one should be 
careful with the singularities at the positions of Bethe roots. This gives 



This equation is also equivalent to the corresponding combination of TBA equations for fermionic 
nodes as it is shown in appendix C.2. Note however that, as we see from (3.1), the functions Yi^i and 
1/^2,2 are not independent but rather they are the values of the same function on two consecutive 
Riemann sheets. So the ratio or the product of (5.24) and (5.30) are enough to fix any of them (up 
to a sign). 

Equation for U: One can find U from the left equality in (5.27). Again we get a finite difference 
equation which can be solved up to a constant A as follows: 

log[/ = logA + log-^ + ^1'*Pc- (5.31) 

Similarly to (5.20), (5.31) is not sensitive to the choice between t^c = log ^1^!*!^°) and 
Pc — log C(w + iO)C(M — iO) which differ by a function holomorphic in the lower half-plane. Nu- 
merically, it is more convenient to use Pc in (5.31) because it quickly decreases to zero when u — > oo. 
We fix the constant A in (E.21). The yet unknown function h is determined below. 

5.2.3 Equations for p and p2 

From (3.13) and (3.16) we get an equation expressing the resolvent G, or equivalently, the density p 
through the Y- functions Y\\ and 
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i + i/i"2,2 _ (1+^f ?p-ip)(i + ^r 



i + -(i + ^+sp+ip)(i + ^-sp+ip) ' (5-32) 

where * p denotes the principal part of the convolution along the interval [—2(7, 2g\. Expressing p 
from the terms without convolution in the r.h.s. of (5.32) we can easily determine p numerically, by 
iterations, as a function of ^^^1^'^-'^ 



To fix Pi we use a similar ratio of Y- functions which can be expressed only through the T-functions 
of the upper band, namely 

1 + ^2,2 ■'^ii-'^iZi ^g_33^ 



1 + 1/^1,1 ■n,-!-^^^-^^:^ 

where each T-function can be expressed through q-functions according to (5.6)'^^. The explicit expres- 
sion for the r.h.s. of (5.33) is similar to the r.h.s. of (5.32), up to the substitution p — > p2 and to a 
number of additional terms involving and 54. As a consequence, p^ is expressed as a function of 
and q3,q4. Note that and (74 are relatively small and vanish in the large volume - the fact 
which is important for our numerical iterative procedure described below. 



^For instance, one gets ^+ = 1+^+ * (p.-gl+Ol -^0] _^[+0] _[-0, 1 ^^^(^[+2, _ 1 ^[+0] ) _^0, ^(^[+2, _ 1 ^[+0, ) _^0] ^ 
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5.2.4 Equation for h 

To obtain this last equation needed to complete our FiNLIE, we employ a couple of Hirota equations 
in the T-gauge: 

T+2T2"_2 = Tr3,2Tl,2 + T2,lT2^3 , (5.34) 

T+iTr,i = Ti,oTi.2 + T2,iTo,i . (5.35) 

We notice that in this gauge To,i = 1, f2^s = t^+"lt^7/l, ¥3,2 = -^"+¥2,3, Ti,o ^ -^[,1^^2,2.7^+ and 
the magic and mirror functions are related for u £ Zq hy ¥2,3 — T2,3, ir2,2(M ^ ^) — '^2,2{u ± f T *0), 
Ti_i(w ± |) 1ri,i(u ± I T iO). Using all these properties and excluding T2,i from both equations we 
get 

i^7/"°lf^+^+°l + -7^+Ti,2 = f [-i+olfl+i-o] + ri,ir2.2-7^+Ti,2 , ueZo. (5.36) 

Inserting the explicit parameterization (3.17) Ti^s = /il^^'/i'^'*! (s + G^+^l — G^^"') wc get the following 
equation on h: 

h{u + iO)h{u - iO) ^ ^ (1-^1,1^2,2) ^ ueZo. (5.37) 

P 

Wc know that h has only one Zg-cut. Asymptotically, for large volume behaves at large u 
as u~^~^ (see appendix B), and hence the monodromy of h? should be trivial for the closed paths 
surrounding its single Z-cut. This topological property is unlikely to be changed for a finite size and 
therefore the large u behavior of h should be the same as in the asymptotic limit. Moreover, since in 
the asymptotic limit h has no poles or zeroes we conclude that at least for sufficiently large but finite 
L (or small but finite g) the discontinuity equation (5.37) uniquely fixes it as follows: 

log/i = log x + Z * log — . (5.38) 

2 VP/ 

We remind that the function J- — -y/To^o can be written in terms of Q-functions. Alternatively 
we can use (3.41) to write it in terms of the Y-functions ^1,1, ^2,2 which is a very convenient equation 
for numerics (see appendix E). 

This concludes the derivation of FiNLIE for the exact anomalous dimensions of the sl{2) sector's 
symmetric states in AdS5/CFT4. Our derivation uses only the properties of T- and T-functions 
summarized in section 3.8, which can be either derived from TEA integral equations or postulated, 
from the explicit knowledge of large volume solution given in appendix B (weak coupling solution 
is enough as well), and an assumption, which can be verified in every explicit computation, that 
qualitative structure of poles and zeroes does not change if to compare with large volume expressions. 

6 List of FiNLIEs 

In this section we will collect together the results of previous sections into the full list of FiNLIE. To 
make the formulas a bit more compact we will employ here exponential notations for the convolutions. 
Namely, by definition, for any kernel S and a function / we will define the exponential of convolution 
as 

r^ = exp(S*log/). (6.1) 
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The equation for Yi^i ( obtained as a sum of (5.24) and (5.30)) now reads'^'^ 

The equation for l2,2 is essentially the same since l2,2 is simply the analytic continuation of Yi_i under 
the Z-cut: Y2^2{u + iO) = 1/Yi_i(m — iO). The T-functions in the r.h.s. of (6.2) are given by 

Ti^s = s + JCsip (6.3) 

and the equation for p in terms of Yi.i and l2,2 is 

1 + 1/>^2,2 _ (1 + i^t^p - 5P)(i + ^r*p - hp) 



1 + (l+/^+Jp+ip)(l+/^-ip+ip) 



e[-25,2.g]. (6.4) 



The upper band i^-functions are expressed through the (/-functions and the LR wing exchange "gauge" 
function U{u) using (5.6) and (5.5): 

=^a,+i — 9l 92 + 92 9i + 93 94 + 94 93 7 lo.oj 
C7 _ „[+a]-[-a] [+a] [-a] [+a] [-a] [+a] [-a] [+a] [-a] [+a] [-a] „x 
=^Q,0 - 9l2 9l2 + 934 934 - 9l4 9l4 - 923 923 ~ 9l3 924 ^ 924 9l3 ' (O.bj 

=^a,-i = (c/[+"li/[-''l) , (6.7) 

with all the g-functions being parameterized through a base of 5 of them: 

91 = 1, q2=P + IC*p2 , 912 =Q, 9123 = C/^, 9124 = ^/^92, (6.8) 

where for the Konishi state we introduce two polynomials Q — {u — ui){u + ui), P = —iu.^^ The 
value of ui is determined by the condition that =3^"^ should have zeroes at position of the Bcthe roots. 
The other q's can be found through the set of Pliicker relations [20, 36] 

909»j = (itoj - 9/9r, (6-9) 
qijkqi = qtjqik - qth^hy (6-10) 

The density p2 in (6.8) is represented as P2 ~ P2 ~ 93^°'94 ^ 94^"'93 where p2 has a finite 
support on Zq and the remainder is a small correction fixed from the consistency with (6.9) and (6.10). 
P2 is then expressed through li^i and 12,2 as follows: 

1±I^ = -^^f^^'^'^^-" , u e [-25, 2g] , (6.11) 

i + i/ri,i ^1+1.^171.^3.2 ^ y' ^ I 

where 2^a,2 = (70+°' 4"°'. 

The function U{u) is defined through (5.31) with the spectral densities expressed through Y- and 
T-functions by means of (5.13), (5. 20), (5.38): 

2;^-2r^^^r2^2.^i,o V / V'^i:i'^i,2i'2,2/ 



^^The overall minus sign in (6.2) is not visible from Vi 112.2 and Vi, 1/12,2- It is chosen so that li^i and 12.2 are 
positive on [—2g, 2g\. The square root ^ Ag'^ — v'^ in Z{u, v)* is evaluated slightly above real axis. 

For the case of more than two symmetric roots we expect these polynomials to be Q{u) = njiiC*" — '"j), ^(") = 
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Figure 8: Algorithm for our FiNLIE implementation. 



where A is defined by relations (E.20),(E.21). 

The M Bethe roots are fixed by the exact Bethe equations which can be written as: 

In (6.12) and (6.13), T is obtained from (E.16) as 

M 

= (yia(«)r2,2(«)) * '°"'''<"-.:'+=-^"<"' J|sinh(7r(u - u,)) , (6.14) 

i=l 

where Ajr is just a constant expressed as (E.18) in appendix E.4. 

Finally, the energy of the state can be then found from the large u asymptotics of the product of 
fermionic Y-functions: 

log Yi,ir2.2 = — + 0(1/m') . (6.15) 
u 

7 Numerical implementation of FiNLIEs 

Using the system of FiNLIE derived in the last two sections, the spectral problem can be solved 
iteratively for the densities of the parameterization (3.16,5.1). We performed the numerical computa- 
tions for the Konishi operator characterized by two symmetric Bethe roots u\ — — it2. We expect the 
generalization to the other operators, with a larger number of roots, to be straightforward. 
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Let us denote a,s X = {p, U, W, {uj}, {%}) the set of parameters characterizing a solution which 
we want to find numerically. It consists of two densities on a short Z-cnt (the density p defined in (3.9) 
and the density p2 which parameterizes W as in (5.12)), one function U (which can be parameterized 
from knowing only its value on the real axis^^), a set of Bethe roots {uj}, and another set of roots 
{ui\ of the polynomial Q from (5.1) ( for the Konishi state i = 1 or 2, ui = M2 and ui = —U2). 

The numerical solution of the FiNLIE system is achieved using the fixed point approach: first, 
the equations are put into a form X = F{X) which allows us to build a sequence of approximations 
Xn+i = F(Xn) converging to the exact solution. The solution is then reached by a large number of 
iterations, each of them consisting of two main steps (see figure 8): 

• On the first step, the T-functions and the fermionic Y-functions Yn, Y22 are computed for given 
densities and Bethe roots: 

1. The upper-band Q-functions are expressed in terms of U, W and Q using the Pliicker 
relations (6.9,6.10). It is clear indeed that (6.9,6.10) allow us to write all the Q-functions 
in terms of the five functions given in (6.8), (by essentially the same steps as at the end of 
section 4.3). Numerically, the main tricky point is to invert the difference equations: for 
instance (713 is obtained from 

u'^=(^-^y -(2rY , (7.1) 

9l29l2 V912/ \Q12J 

/ 2 \ [2fe-|-l] 

The l.h.s. behaves as u^L--f-2M .^^rhen u > 1, hence = V, ( +' _ ) , where 

912 ^ ^k>0 y q+^q^^ j 

the sum converges very fast. It is also possible to rewrite this equation as an integral 
equation, similarly to (5.19). At the end of this first step, once all these Q-functions are 
obtained, one computes =^'s using the equations (6.5-6.7); 

2. The right band T-functions are expressed in terms of p, using (6.3); 

3. Then it is possible to compute Y\.\ and 12,2 using equations (5.24,5.30). To do this, we use 
the 5^-functions, the T-functions, and the position of Bethe roots. 

• On the second step, the set of parameters X is expressed from these T-functions and the fermionic 
Y-functions: 

4. First, Uj is obtained by the requirement that S^\_^(uj + j/2) = 0; 

5. Then the density p is extracted from Y\\ and I22 using (6.4). In the same way W is 
extracted from Yn and F22 using (6.11); 

6. The function U is then found from (6.12), with T fixed from (6.14); 

7. Eq. (6.13) is used to express the positions of the Bethe roots {uj}. 

Finally, for the initialization point Xq one can either use the asymptotic solution, which works 
perfectly for small (7's, or extrapolate X from smaller values of g. 

Numerical results: Our, rather preliminary, numerical realization of FiNLIE's gives for the Konishi 
state the results written in table 3. 

^^The large-n behavior of U is given by C/ ~ -^j. — ^—i)!'^ (see (B.18) and the discussion at the end of appendix B), 
allows us to write a Cauchy-kernel representation for it as follows 

1/2 = -L_^ — ^-^Av , where p„ (u) = 2 Re f C/^ ■ x^+^-l) 
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9 




Et^a [22] 


-E^FiNLIE 


numerical error 


0.25 


4.6154 


4.6147 


4.615(7) 


10-3 


0.5 


5.7362 


5.7120 


5.712(9) 


10-3 


1 


7.7285 


7.6044 


7.60(53) 


10-3 


1.6 


9.5749 


9.3874 


9.38(23) 


5 • 10-3 



Table 3: Numerical energy of the Konishi state, for a few values of the coupling constant g. The 
outcome of our FiNLIE iterations is compared to the previous TBA-iterations [22] and to the prediction 
of the Asymptotic Bethe Ansatz (ABA). 

These results are in perfect agreement with earlier numerical study, using TBA approach, of the 
Konishi state [22] (confirmed by a similar numerical study of TBA in [25]): when g < 1, the numerical 
precision is essentially the same as in [22] . 

For a given state which we take as Konishi state here, our numerics showed the convergence of the 
algorithm3^ to the exact solution of all the above-written equations, and this solution coincides with 
the previous numerical studies from the original TBA equations. We notice however that numerically 
it is more efficient to use the old exact Bethe equation derived in [7] rather then (6.13). 

In figure 9 we can see that the Y-functions Yafi, Yi i and Y2.2 obtained from FiNLIE iterations 
(dots) coincide with a good precision with the data of TBA-iterations (solid curve) . We can also see 
a sensible deviation from the asymptotic limit of these Y-functions (dashed curve) . 

For Konishi state, the plots of various densities resulting from the iteration procedure are presented 
in figure 10. The densities p, p2 and pu (black curve) are compared to their asymptotic values 
(B.5), (B.13), (B.20) (dashed gray curve). We see that at 5 = 1.6, the densities are already quite 
different from their asymptotic value. In particular, figure 10b shows that the density p2 sensibly 
deviates from its asymptotic semi-circle shape. 

In conclusion, our preliminary numerical results confirm the validity of our FiNLIE system. We 
plan to improve the precision and speed of our numerical procedure in the near future. 

8 Conclusions 

In this paper we have shown that using the integrable Hirota dynamics, together with a few relatively 
simple and natural assumptions about the analyticity properties of T-functions (parameterizing the Y- 
functions) we can transform the AdS5/CFT4 Y-system into a finite set of non-linear integral equations 
(FiNLIE). 

Our FiNLIE remotely resembles, at least in spirit, the famous Destri-de Vega equations known in 
the literature for some other 2D QFT's. But the method of derivation, first proposed in the context 
of the SU{N) principal chiral field (PCF) in [4], is conceptually very much different from the original 
approach of Destri and de Vega3''. The relation between two approaches still needs to be clarified, 
though some interesting results on this way are obtained for the principal chiral field model in [43] . It is 
relatively easy to formally solve the integrable Hirota equations in T-hook in terms of a few "boundary" 
functions, which is done by means of the Wronskian determinant formulas [18, 35]. However, it is much 

^®For Konishi state, it is checked numerically at least in the range < 9 < 2. 

^^though the resulting equations can be identical in some simplest cases, like Sine-Gordon or Gross-Neveu models, see 
[4]. 
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(a) Middle nodes Y-functions (Va.o) (b) Fermionic Y-functions {Yi i and I/I2.2) 



Figure 9: Numerical Y-functions for Konishi state at g=1.6 : The Y-functions obtained from 
FiNLIE iterations (dots) are compared to the outcome of TBA-iterations [22] (sohd hnes) and to the 
ABA expression (dashed Hnes). The figure 9a shows Yi_o (black), ^2,0 (blue), Y^^ (green), ¥4^ (red) 
and Y^fi (violet), while the figure 9b shows Yi^i (green) and its continuation 1/^2^2 (blue). 




(a) Density p when g = 1.6 (b) Minus density p2 when g = 1.6 (c) Density p^j when g = 1.6 



Figure 10: Numerical densities obtained for Konishi state by our FiNLIE algorithm: 

In figure (a) the density p describing the right band is plotted (black curve) and compared to its 
asymptotic expression (B.5, dashed gray curve). The same comparison is presented in figure (b) for 
—p2 and in figure (c) for pu. 

more difficult to establish the right analyticity properties of those functions: asymptotic conditions, 
analyticity strips, analytic continuation through Z-cuts present in all T- functions, etc., and thus fix 
the full physical solution. 

As for the case of the PCF, we observe that the analytic properties of T-functions can be reasonably 
simple (i.e. exhibiting well established analyticity strips) only in a well chosen gauge which is different 
for right, left and upper bands of the T-hook. 

Our fundamental observation, mentioned already in [18] but clearly understood and employed only 
in the current paper, is the fact that the Z4 symmetry of the string coset model, explicit in the classical 
system, can be promoted to the quantum level as a symmetry under the analytic continuation of the 
T-functions w.r.t. their representational indices (a, s). For a given T-function, this is only possible 
on a certain sheet of the Riemann surface having only short Z-cuts. We call it the magic sheet. 
The Z4 symmetry, the quantum analogue of unimodularity (quantum determinant = 1), and a few 
other natural assumptions allow us to fix the analytic properties of each T-function in the appropriate 
gauge, as well as the transitional gauge functions relating T's for three different bands of the T-hook. 
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Only then, knowing the whole structure of all T- functions, we have enough of analyticity input not 
only to parameterize all T-functions and the gauge transition functions but also to constrain them by 
additional Riemann-Hilbert-type equations, thus closing the whole system of FiNLIE. 

Our system of FiNLIE (summarized in section 6) shows that we achieved our conceptual goal of 
getting a finite system of equations out of a set of analyticity properties of T-functions. Probably 
the actual realization is still perfectible and we only begin to see the hidden "simplicity" of the whole 
problem. It is even conceivable that there exists a "quantum spectral curve" of an infinite genus, with 
infinite number of sheets connected by Z-cuts, uniformizing the analytic structure of all basic functions 
(presumably Q-functions) of the model at once. 

We also demonstrated here that our FiNLIE is a useful tool for the computation of the spec- 
trum. Our first, preliminary numerical implementation of the FiNLIE (which, with a certain effort, 
can be certainly improved in efficiency and precision) allows to check that, within reasonable error 
margins, our numerics reproduces the known results [22], thus perfectly confirming the correctness of 
our FiNLIE. We also proved here the equivalence of our FiNLIE to the TBA equations of [14-16, 22] 
analytically (the paper [11] was especially useful for that). 

Our method rendering a finite system of equations for the AdS/CFT spectrum can be certainly 
generalized to any state of the model but the details for the other states still have to be worked out. 
We also hope that our FiNLIE will allow in the nearest future an efficient way for constructing the 
systematic weak coupling expansion. It will be more difficult to do the same for the strong coupling, 
thus making it possible an efficient higher loop calculation in the string sigma model, but our FiNLIE 
could hopefully provide us with some clues also for that. It is especially important in view of recent 
analytic results obtained by a very different method at strong coupling [44-48]. Indeed FiNLIE, unlike 
the TBA, provides an extensive knowledge about the analytic properties of underlying functions on 
the whole magic sheet. It should be also easier now to attack another important limit of the theory 
— the BFKL approximation. Hopefully FiNLIE will allow to re-derive the known leading [49] and 
next-to-leading [50] BFKL approximation for N=4 SYM and to attempt a systematic study of this 
expansion in N=4 SYM. Another important quantity which would be interesting to compute is the 
"slope function", the exact form of which was conjectured recently in [51]. 

Finally, we hope that our result will bring us closer to the understanding of the mystery of the 
AdS/CFT correspondence. At least in technical terms, due to FiNLIE the gap between two sides of 
duality seems to be much narrower now than before. But an important ingredient - a derivation of 
the Y-system and of FiNLIE from the first principles, i.e. directly from both the string sigma-model 
and the N=4 SYM theory - is still missing. 
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A Z4 symmetry and strong coupling limit 



In this appendix we remind the classical Z4 symmetry of the string coset at strong coupling and 
motivate its quantum generalization given in the main text. 

We know from the group theoretical analysis of finite gap solutions for the superstring on 
Sp^^)x'sp(4) ^^^^^ [1^! 52] that in the classical limit our T-functions become characters, explicitly 
written in [18] (see eqs. (4.12-21) there), in the highest weight unitary irreps a" of PSU{2, 2)4): 

Ta,s = tTa^si^iu) , (A.l) 

where fl(u) is the classical monodromy matrix with the eigenvalues (/ii, /i2, fJ-s, fJ-il^i, A2, A3, A4). Note 
that this identification was made in [18] in the mirror sheet and the expression (A.l) should be con- 
sidered in the mirror kinematics. The classical monodromy matrix as a function of spectral parameter 
u has only one large cut {—oo,—2g] U [+2g, +00) with an essential singularity at the branch points 
u = ±2g. 

The Z4-symmetry of the coset sigma model [24] imposes the following relations among the eigen- 
values [23]: 

Mi(7/) = 1/M2(«^), M3(«) = 1^4(w^), Xi{u) = 1/X2{u'^), X3{u) = 1/X4{u^) . (A.2) 

Here we denote by f{u^) the result of the analytic continuation of a function / following the full circle 
around the branch point u ~ 2g along a path avoiding square root cuts of a finite gap solution which 
are present in the eigenvalues but are absent in the monodromy matrix'^^. 

For classical T-functions, or characters, given by the formula (2.19) of [17] the property (A.2) 
implies the following symmetry 

TaAu) = {-lyTa.-siW), if |s| > a 

where the functions Ta^-s represent the analytic continuations of the functions Ta,s with respect to 
the argument s from the values s > a. _s should not be confused with Ta.-s entering the Hirota 
equation (1.1) on the T-hook. r_a,s('i*) is also simply the analytic continuation oiTaA"^) from a > \s\ 
to negative a. 

The classical limit has a rather degenerate analytic structure w.r.t. the full quantum case, with the 
eigenvalues (and hence the T-functions) having only two branch points at ±2f/. In the full quantum 
case, we have, a priori, an infinite system of branch points, at u = ±.2g ± i{ -|- n), n = 

0,1,2,3,.... At strong coupling — >■ 00 there is no way to distinguish between the branch points 
with different n's. Hence, the quantum version of (A. 3) is potentially ambiguous. However, discussed 
in this paper analytic properties of Q-functions give a natural suggestion how the proper quantum 
version should look like. 

Let us consider the right-band T-functions and use the Wronskian parameterization for them, e.g. 

\+s\ [—si [+sl f— s] 

Ti.s = qi P2 I2 Pi ■ The transformation Ta.s —> Ta,-s simply becomes the shift of the spectral 
parameter'^^ u in the mirror kinematics: gj"''*' ~^ ?! p\ ~^ p!^"^- N^ow we consider only g-s and 
note that combination of the shift with analytical continuation clockwise is nothing but the shift 
qj^'*' avoiding the short cut [—2g, 2g\, see figure 11. The same observation holds for p-s if the 



^^these cuts may be interpreted as a condensate of Bethe roots. 

^^Note that for (/-functions sucli a shift is not neghgible, even in the character hmit. 
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Figure 11: Equivalent representations of transformation of g-functions following from 
(A.3) 



continuation u"' is made counterclockwise'^^. Therefore (A.3) for the horizontal band is reformulated 
as: 

TaAu) = i-irfa,^s{u), (A.4) 

where Ta^-s is defined by the Wronskian ansatz for the right band considered at —s and with all shifts 
of q- and p- functions avoiding the short cut. 

Generically, Q-functions on the quantum level have infinite number of cuts, which makes the above 
arguments very ambiguous. However, in the T-gauge we know that the Q-functions can be chosen 
to have only one short cut. This makes the T-gauge (and any other gauge with only two short cuts) 
more suitable for the generalization of (A.4) to the quantum case. That is probably why one gets 
(3.46c). A similar analysis for the upper band is more tricky. However, in appendix C.4 we show that 
upper-band Z4 symmetry can be derived from the right-band Z4 symmetry. 

The outcome of our observations is that the quantum Z4 symmetry (3.44c) and (3.46c) is the 
most natural generalization of the classical Z4 symmetry of the string sigma model on so{i^5)>fs)l6) 
coset. However, we should note that it still remains to be proven that the generalization we propose 
actually reduces to the classical Z4 symmetry in the strong coupling limit. One should carefully 
analyze the possibility of interchanging the order of making analytical continuation and taking the 
strong coupling limit. Also note that in the quantum version we use two different gauges for the upper 
and the right band to formulate the Z4 invariance, whereas in the classical case they should coincide 
up to a constant for |Re(u)| < 2g. This effect is accounted by the i-periodic function T relating two 
gauges which in the scaling regime used for the classical approximation is indeed approximated by a 
constant for |Re(u)| < 217. 



B Large volume asymptotics 

The asymptotic large L solution of the Y-system is used in this work as a starting point for the iterative 
numerical solution of FiNLIE. The positions of cuts/poles/zeroes and the large u expansion for the 
asymptotic solution are explicitly known. We use this information to deduce the qualitative structure 
of cuts/poles/zeroes and to derive the large u expansion of the exact solution. 

We consider only the s[(2) sector, hence all the Bethe roots are expected to be real. We also 
study here for simplicity only symmetric configurations of these roots, i.e Uk — —UM~k with an even 

There is no apparent contradiction in using the opposite directions of 7 for continuation of q and p and the same 
direction for continuation in (A.3). The reason is that 7 defines a monodromy of order 2 as it follows from (A. 2). 
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number of roots M . Some of the formulae written below are simplified by using these assumptions. In 
particular, the vanishing of the total momentum is already satisfied due to the relation*^ J^j, ^ = 1. 
It is useful to define the following functions 

P = ^^^1 + , S = - . (B.l) 

Note that P — + ... is a polynomial in u, and is a polynomial times pure square root. 
B.l Right band 

An expression for the large volume limit of T-functions was presented in [7] . On can check that up to 
a gauge transformation this expression for the case of the right band is equivalent to: 



/ 11 + dl+^'l 11 + $["''1 \ 7a, 



where is the asymptotic value of the anomalous dimension and $ = (-l)^^/^fpy is defined as a 
function with short Z-cut. The normalization constant Cj- is chosen so that 7i,s — > s, u — c». For 
configurations of the Bethe roots without poles in the analyticity strip of 7i,s one can write a spectral 
representation for 7i,s: 

fi,s^s + IC,ip. (B.3) 

The density p is given by: 

$[+0] _ $[-0] 25 

P = Cr = Cr—, , ueZo- (B.4) 

(l-$[+"l)(l-$[-ol) ^Q+ + Q--2P' ° ^ ' 

The above expressions give indeed the large volume limit in the T-gauge introduced in section 3. One 
can see that the denominator in the rightmost expression of (B.4) is a polynomial of degree M — 2. 
For the case of two symmetric magnons (M = 2) (B.4) simplifies to 

p = i^/MEZ = - - (B.5) 

7as +4 xT - X. 



and (B.2) simplifies to 



where we choose ui = —U2 > 0. 



ii,s = s -T , (B.6) 



B.2 Upper band 

For two symmetric magnons there is an expression for the T-functions in the =f7-gauge which is very 
similar to (B.6): 

1/xI'^l - l/x[~"l 
l/x+ - 1/x^ 

'^^In the generic case the momentum is given by (2.9). 
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In this form it is easy to see that ^i^i{uj) = 0. For general M the asymptotic expressions for the 
c^^'-gauge can be read off from the Wronskian parameterization in [18]. The explicit asymptotic values 
of our basis of Q-functions look as follows: 

91 -1 , W = q2 = C.^(M-S^ , qi2 = Q = Q- (B.8) 

Here Cy = - — M is a polynomial solution of the following equation: 

M+ - M- = 20 - P+ - . (B.9) 

The additive constant cannot be fixed from (B.9) but it is irrelevant due to the H-symmetry (4.3). A 
convenient choice for this constant is given below. 

In the leading large volume approximation the T-functions in the ^-gauge are given by 

= I4^[+-l - , #,,0 = O'+^lQ^-'^l , = qI+JqM]' ■ (B-10) 

One can see that 

AAu,) = 0, (B.ll) 

so the poles in the denominator of S^a,2 cancel. The normalization constant introduced in (B.8) 
is chosen so that = au^'^~'^ + . . . for large u. 

M — 1 

Since W has no poles and for large u it behaves as one can construct the following spectral 

representation: 

W = Pm-1 +IC*P2 , P2 = -2C^S , (B.12) 

where the zeros of the polynomial Pm-i = ~ m- 1 " H jlT ^ ~ ^i) fixed by condition (B.ll). Note 
that we used a freedom in the choice of the additive constant of M to constrain the form of Pm-i- 
For AI — 2 we get simply 



1/x X ~ 1/x 

l/x+ - 1/x^ ' " l/x+ - 



W = -iu- ^ ^ , _ , P2= , , + ■ (B-13) 



Expression for U: Comparing 



^2 /jj[+a]^l-a]\2 

Yafi = ^^ . , lmu>a/2 (B.14) 



to the asymptotic expression for Y^.o in [7] we find 



U{u) ^ ^=B^-\u)^)^ , lmu>0, (B.15) 

U{u) = ^i?(+)(u)— y , lmu<0, (B.16) 
VC> V{u) 



where Au is a real number and y is a solution to the equation 



jj_ 0-BEs(u,U/c) = „/_■)_ _[]_ Crniirr(u,Ufe) . (B.17) 



V- ~ Bi-)+ -IJ- - I B(-) 

fe=i \ fe=i 
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The function fimirr first time appeared in [15, 53] and was later identified in [34, 54] with solution of 
the mirror crossing equation. It has a particularly nice behavior under fusion: (t(u, v)^°'^° as a function 
of u has only two cuts, Za and Z-a, on the mirror sheet. As a consequence, V has only one cut, Zq, on 
the mirror sheet. From integral representations for a^irr given in [15] or [34] one derives the following 
explicit expression for V: 




Xo[u,v) = log a;, 

y 

where x = x{u), y = x{v), Xu: = xiuj), Xq = x{q). 

Note that log 1^ is a purely imaginary function in the mirror kinematics (which is consistent with 
(B.15)). The large u behavior of V is defined solely by xo- One gets logF ~ — ^ logu which implies 

\ogU -^^^logu , u->oo + iO. (B.18) 

Asymptotics of h: The large volume expression for h is rather complicated. However, we would 
need only the large u asymptotics of this expression. From the condition Ta,2 = "^2^ one can relate 
h to the other, already known functions: 

where C is some real number. One can see from (B.23) and (5.18) that / ~ u^^^ and thus from the 
previous relation one has 



L + 2 



u~^-^ (B.20) 



or 

h^ir^i" . (B.21) 



We expect that in a finite volume this asymptotic behavior remains unchanged, see discussion after 
(5.37). 

Some useful expressions: It is known that the product of fermionic Y-functions is very simple 

^^^^^^^^ - BWbT) ' (^-22) 
so that from (5.13) one can easily determine the large volume expression for B : 



B=(^) . (B.23) 
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The function T can be determined by its discontinuities and zeroes. By applying arguments of ap- 
pendix E.3 to the large volume expression (B.22) one gets: 



T^-K fr '''''' 

11 [2a-l] [2a-l] ' (^-24) 
a=-oo ^ 1,1 -"^2,2 

The exponential numerical factor in the numerator is needed to ensure the convergence and A^? is a 
normalization constant. 

Large u behavior of exact quantities: Above we derived the large u behavior of different quanti- 
ties in the large volume approximation. The result however remains the same at finite volume provided 
that we use the exact value for the energy (2.9) and anomalous dimension. This was already shown for 
^1.1^2,2 (see Eq. (3.40)) and for h. Finite volume corrections also do not change the leading large u 
term of ^,s>o and Ta,s- Indeed, the magnitude of corrections in these quantities is defined by If^lJ^" 
for some positive integer a which is smaller than (il+^l^I""!)^-^/^. From the definition (5.13) for B 
and / we can now conclude that / ~ 

and B ^ 1 — ^ at any volume. Finally, we can invert 
relation (B.20) to prove that (B.18) holds at any volume. 



C Equivalence of FiNLIE to the TBA equations 

The finite volume AdS/CFT spectrum problem was previously analyzed in the literature using an 
infinite set of the TBA equations. These equations were proposed in [14-16], following the discovery 
of the Y-system [7], and passed a number of important checks [17, 22, 25, 44-46, 52, 55-57]. It 
is important to verify that FiNLIE is equivalent to the TBA equations and this is the goal of this 
appendix. The TBA equations can be also viewed as a departing point to derive the properties of the 
transfer matrices which we summarized in section 3.8. Alternatively, these properties were considered 
there as the initial assumptions for the derivation of FiNLIE. 

The TBA equations and the notations we use are the ones from [15, 22]. The papers [14, 16] 
contain an equivalent set of TBA equations in different notations, but without the driving terms for 
the excited states. 

C.l TBA for Yi,s and ^^,1 

We start our analysis by the TBA equations for the right band: 

°° 1 + Y 

iog(ri,„+i) - V/c„,™*iog(i + i/yi,„+i) + /c„siog— — ^ , n>i, (c.i) 

ICn,m = - D-2) [n]u[m]ulC , D = e^^" . (C.2) 

The fusion operator [s]d is defined in section 2.2. 

Now one replaces the Y-functions with a combination of the T-functions using 1 + 1/11. m = 
7-+ T — 

7j= — — . This leads to the chain cancellations of an infinite number of terms in the sums leaving 

+ l /l,m-l 

US only with a finite number of terms in the r.h.s of (C.l). We call this procedure "telescoping"^^. 



^The same kind of simplification was also used in [12, 13] under the name of "chaining relations" 
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Indeed, after the substitution of T-functions and the shifts of the contours of integration in the terms 
containing logTj^^^ one gets: 

ICn^m * log(l + l/Yi^rri+l) = 

/ oo oo oo \ 

(D+2-D-2)[„]d K](D + D-i)[m]D " ^[TO+l]D-5I["^-l]D U:*logri^™. (C.3) 



m=0 



(C.4) 



Since (D + D ^)[to]d — [m+ 1]d — [to — 1]d — 0, aU the terms in the infinite sums above cancel out 
except the boundary ones. We therefore get: 

oo 

^ K-njn * log(l + l/Yl,„i+l) = log(Tl,„+2Tl,„) - (/C„+i + K-n-l) * log(Tl4) 



m— 1 



iog(ri,„+2ri,„) - /c„ * iog(ri'+'"°'ri'T'+°i)- (c.5) 



The formal manipulations that we made are valid only if the integration and summation can be 
interchanged and only if the contours of integration do not hit poles or branch points when being 
shifted. The T-gauge given by (6.3) enjoys these nice properties and that is why we used it for the 
derivation of (C.5). Note that the r.h.s of (C.5) is not invariant under an arbitrary gauge transformation 
while the l.h.s is; that is so because the telescoping procedure can be performed only in a specially 
chosen gauge. 

One can substitute Yi_n+i = ^r^^+i" ^'^-^ TEA equation (C.l). The numerator is 

cancelled against log(7i,„+27i.ri) term in the telescoped expression (C.5). Therefore this TEA equation 
is reduced to 

- iog(r2^„+i) = -/c„ * iog{Tl^,'-°^T};,'''°^) + /c„ s log ^ + 



1 + 1/Y22 

/c„ * iog(ri'r""'r/r+'^0 + /c„ i log 



-7— q-+ 
'2,2 '2,2 



^2,3 



= -/c„* iog(ri;iri.i) + /c„ * log ■ (c.e) 

'2,2 '2,2 

Using the representation (3.10) for 72,n+i we get: 

l0g(r2,„+l) = l0g(l + G["+21 - G["1) + l0g(l + G[^"-2] _ Ql-n]-^ 

= ^„ * log(l + G^l - G[+ol)(l + Gl-21 _ Gl-ol) = /C„ * log J^-^ . (C.7) 



72 



2,3 



T 1 — 

Note that the expression in the argument of the log is equal to 72, 1 for u > 2g and ''.^ '^•^ for u < 2g 
so that (C.7) reduces to: 

^n*log;^^^ =0. (C.8) 

This equation should be valid everywhere on the real axis and for any positive integer n, which is only 
possible if 

%\f,:,^f2,i. (C.9) 
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From the Hirota equation for 7i,i we have T^iTi i ~ T2.1 + 7i,o7i,2, and consequently*'^ 

ri,o = gi+°i + g[-"i = 0. (C.IO) 

We thus derived from the TBA equations the statement (3.15). 

The derivation presented here can be reversed. Hence the set of TBA equations for yi,s>2 is equiv- 
alent to the condition 7i,o = 0; with an additional assumption that 7i,s S As and with a large volume 
asymptotics Ti.s ^ s + 0{l/u). The condition of the absence of poles and zeroes in the "telescoping 
strip"** — 1/2 < Im(u) < 1/2 is also needed, however it is relevant only because we consider here the 
sl(2) sector. In case of a general state zeroes may be present in these T-functions and they should 
lead to the appearance of additional driving terms in the TBA equation (C.l). 



Equations for "pyramids": Let us now consider the TBA equations for {a, 1} nodes of the T-hook: 

log(y„+ia) = [n-l]Dlog- 

— ICr. 



Bi+) Ri+} 
n+l]Dlog 



* lofffl + 




Mn+l,m * l0g(l + Ymfi) 
1+>1,1^ 



n> 1. 



(C.ll) 



V2,; 



Here and below the summation sign over the doubly repeated m from 1 to 00 is systematically omitted. 

We want to perform the telescoping procedure similar to the one for the right band. For that it is 
natural to use the ^-gauge which has the simplest possible large u behavior. One should be especially 
careful with the terms containing Ym,o in the r.h.s. which have a complicated kernel*^. Moreover the 
term with fermions does not have a form, which is easily expressible in terms of the T-functions of 
the upper band. Both complications can be overcome at once by subtracting from (C.ll) the TBA 
equation for Yi i after applying to it /C„(S): 



/C„® (^logYi,! = log 

The result of subtraction reads 
n[+n+il 



7e("™) * log(l + Y^,o) ~ ICr,. * log 

-L + -1-/ il,m+l 



logi; 



n+1,1 



log 



Q[-n+l] 



JC 



.1 l0g(l -I- Ymfi) - ICn,m. * l0g(l -I- Ym+1,1) + /C„ S log 

(D -D-i)[n]DHD/C . 



1 + >2,2 



where ICf^,„,_i 



(C.12) 



Now one can perform the same steps as for the right band. The difference is that ^a.s functions 
do contain zeroes inside the telescoping strip. However, the only fate of the residues coming from 
these zeroes is to cancel the driving term log q[„-i ] ■ Note that the driving term is definitely such only 
for Konishi-like states at a sufficiently small coupling. In general this term may be more complicated 
[22, 58]. In any case this term is constructed precisely in the way to cancel residues coming from the 
contour shifts. Therefore the result of telescoping procedure is not sensible to the structure of the 



*^The property (C.IO) can be equivalently reformulated in terms of Y-functions as ^7*2^l~2 = + ^1,3)- The latter 
equation was first time derived from the TBA system in [16], however the magic properties of the T-gauge following 
from (C.IO) were not discovered there. 

**This is the strip in which the contours are shifted. 

*^we use the notations from [7] 
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driving term and in this sense it is more universal than the TBA equations. In complete analogy with 
the right band, the final equation we get for the upper band is: 



= 0. (C.13) 

Again, the equivalence works in both directions: one can derive the set of TBA equations for Ya,!, 
a > 1 from (C.13). The required assumption, in addition to (C.13), is the existence of the ^-gauge 
with S^a.s G -4q-|s|+i ^iid a polynomial behavior at infinity. 

C.2 TBA for Fi i and ^2,2 

One can make the telescoping procedure for the TBA equations for Y\^\^2,2 and Y\^\IY2,2 and get 
precisely the integral equations (5.24) and (5.25). We leave this exercise to a curious reader, while 
here we will use the results of [11] for the derivation of equivalence between (5.24), (5.25) and the 
TBA equations. In [11] it was shown that the fermionic vacuum TBA equations are a consequences 
of the discontinuity conditions: TBA for Y\.\Y2.2 is equivalent to the discontinuity condition (1.7) 
of [11] and the TBA for Yi, 1/1^2, 2 is equivalent to the discontinuity condition (F.5) in that paper^^. 
The discontinuity conditions are proved to be correct also when the excited states are considered [13]. 
We will show now that these discontinuity conditions are equivalent to the analyticity of B and C 
functions, in the upper half-plane — the essential property needed for the derivation of (5.24) and 
(5.25), respectively. 

TBA for the product of Y\^\ and 12,2= The condition (1.7) of [11] is written in our notations as 

n 

disc ((iogri,ir2,2)[+'"i) = - ^ disc (iog(i + ri'o"""')) , « > i . (c.i4) 

a=l 

We express l + Fa.o through the T- functions using l + Ya.o = ^ ^i'o5^°'°i • Again, massive cancellations 
occur in the sum and the final result is: 

disc ((iogr,,r2.2)[+^"i) --disc log tL"-,;;] ■ (C-15) 

moreover since S^}^'^^ ^' and q does not have a cut on the real axis we conclude that (C.15) is 

precisely the condition of analyticity of B, defined in (5.13), in the upper half-plane. As we know, the 
analyticity of B allows the construction of the gauge transformation / via (5.20), which then implies 
existence of T as we show below. Let us consider then the LR-symmetric gauge T^^s — \J ■'^a,s-'^a-s 
and perform the following gauge transformation: 

Ta,s - . (C.16) 

Since / by construction is regular above T- functions have the same analyticity strips as 

functions. Moreover from (C.16) one can show that 

Ta^To^ ^ ^ 

2,3-1-0,0 



*®Though (F.5) is claimed in [11] to be a consequence of other discontinuity relations and Y-system equations, we 
found that the derivation of this claim contains a logical gap, at least in the way it is presented. Therefore we consider 
(F.5) as an independent requirement on the discontinuities of the Y-functions. 
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Next, since T-functions are real by construction, the complex conjugate of (C.17) reads J''^ +^ — 1. 

T2,3Tq_o 

This immediately implies T^J^q = Tq"q and T3 2 — T2,3, and thus the T-gauge is indeed the gauge 
introduced in section 3 (Z4-symmetry in this gauge is shown below) . We therefore arrive at the initial 
assumptions that we used in the main text for derivation of regularity of B in the upper half-plane. 

TBA for the ratio of Yi i and ¥2.2- The TBA equation for the ratio can be reformulated as a 
discontinuity equation (F.5) of [11]. In our notations it reads as 



\ a— 2 ^ / 1 / a— 1 ^ 



= 0. 



After the substitutions 1 + Ya s — ^' "'^ — , 1 + l/li a = t ^'""J^'" — and the telescoping of the 

' a-\-l ,s a — l ,s ' 'l.a + l'l.a — 1 

sums we get 




disc (21og^^-log^-log( ^ ) I =0. (C.l^ 



0,0 



It easy to see that this equation is equivalent to the condition of regularity of C (5.27) in the upper 
half-plane. 

Again, the only assumption we used is the existence of a T-gauge with Ta.s G -^s-a, s > a. Let 
us show now that the T-gauge defined by (3.45) has the same analyticity strips as the i^-gauge. For 
that it is enough to show that the function h relating T- and T- gauges by 

T,,, = (/l[+^l/l[-^l)["l°ra,s (C.19) 

is analytic in the upper half-plane. To show this, we note that from analyticity of C in the upper half- 
plane and (5.27) it follows that is analytic in the upper half-plane. Therefore discontinuities 
6 = A(log ft.!^"!), if non-zero, are the same for arbitrary n > 1. Let us introduce an i-periodic function 
V with A(log'P) = 6. Then h' = h/V is analytic in the upper half-plane. From (/i[+''1/i[-'*1)[21d = 
T2,s/72,s = T2.s/72,s e As-1 we conclude that 2A(log7^7') = and hence we can choose V to be 
such that VP = 1. Therefore h^+'^h^-'^ = h'^+'^h'^-'\ i.e. we can always adjust h, by multiplying by 
irrelevant i-periodic function, to be analytic in the upper half-plane. 

C.3 TBA for Ya,o 

It is again useful to start from the discontinuity condition (1.6) of [11, 13] instead of the equivalent 
TBA equation for Ya.o- This condition reads in our notations as follows: 

disc (disc (log Yi+)P»1) = 2disc( log(l + y[;^i"l) + log(l + if ;-;"!) ) -logYi^ir2,2, 

m— 1 



n > 1. Let us first compute disc (log Yj^Iq) . We use Yi^o = o ~ T^t ^^"^ '^^-^ ^ "^^^^ 

disc (logYj^o) = 2disc (logTj^j^) — 21ogTj'"j^/']rj'"]^, where Ti.i is defined as follows: It coincides with 
Ti^i below i/2 and has a short Zi cut; all other cuts of Ti_i above Zi are of the Z-type (long). 
After telescoping we get 

=- log 1-1,1^2,2. (C.20) 
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Now we use Tq^i/T^^ — J^^/T^-^^. Since the analyticity strips for T-gauge and the periodicity 
of Tq.o = J'^ were already established (from TBA for Yi. 1^2,2) one can use (3.41), which is derived 
indeed from the periodicity of To, and the analyticity of Tio. It states: discJ^[^"+^l = — log Yi. 1^2.2- 
Therefore (C.20) simplifies to 

disc logi'^i+^1 =0, n = T;^. (C.21) 

This means that i has only two short cuts Z±i, and is regular otherwise. This condition is also in 
our list in section_2.2. In the previous subsection we showed that h is analytic in the upper half-plane. 
Since Tij = h^h^ Ti^i we see that h has only one Zp-cut and is regular elsewhere in the complex 
plane. 

Let us now show that h can be chosen to be a real function. For that we use equation (5.37). It 
was derived for real h but this assumption can be relaxed in the derivation: 



/j[+Ol^[-0] ^ 



^+(1-^1,1^2,2) 



P 

The r.h.s of this equation is regular on the real axis. Indeed, p = G^+^l — hence it behaves 

as a pure square root on the real axis. Since yl\°^Yl%°^ = l/(W,'i"'^2,2°') and j:-[i+o1/Jf[i-o] ^ 
l/(^i~i°'^2t°')' numerator of the r.h.s also behaves as a pure square root on the real axis. Therefore 
we conclude that U+^'^h^-^'^ ^ U-^^h^+^y Hence h/h is a function analytic everywhere. Its large u 
behavior should be at most polynomial as it follows from the asymptotic solution, see (B.20). For the 
asymptotic solution one can check the absence of poles/zeroes in h/h and we assume that this is also 
the case for the exact h. Therefore h/h should be a constant. It is always possible to redefine the 
normalization of the T-gauge so that this constant is equal to one. Hence h can be indeed chosen as 
a real function. 

C.4 Deriving Z4 invariance 

Z4 invariance of the T-gauge: In section 3.4 we show that the condition 71 = implies the 
Z4 symmetry of Ta.s and only two magic Z-cuts for 7i,s. As a consequence, 72, ^ is given by 72, s — 
tI\'^^'TI i^"^ and hence it has only four magic Z-cuts. Using the reality of h we see that Since T- and 
T- gauges are related by (3.22), the Z4 symmetry of the T-gauge follows from the Z4 symmetry of the 
T-gauge. Since h has only one cut, the magic cut structure of the T-gauge is also preserved in the 
T-gauge. 

Z4 invariance of the T-gauge: The simplest is to derive Z4 symmetry on the border. Since 
Ta,2 — T2.0 ~ T2,a and T-gauge is Z4-symmetric, we get: 

ta,2 = t_,,2 . (C.22) 

Moreover the TBA equations Yi,^ imply (see (C.13)) 

to,i = . (C.23) 

Then, the followings equalities 

^ = f ^ 1 — = 1 (C.24) 

To,o To, 2 ' T2,o 

and Hirota equation in the magic at the node (0, 1): Ti,iT_i,i + To,oTo,2 = imply that 

Ti,i = -T-1,1. (C.25) 
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Finally, consider u slightly above real axis with Re {u) > 2g, so that f — — . Then: 

/ \ 2 



rp+ rji- \ X-+ TT- / \ 77+ f"- j:+ \ 

-'-o.o-'-o.o \ / y-^ J- -I I 

From the Hirota equation in the magic T(|qTq q = Ti.oT_i,o and (C.26) it foUows that 

ti,o = t_i,o. (C.27) 

Relation between symmetries in the right and upper bands: Let us show that the property 
Tq.i = can be obtained from the property Ti^o = instead of the TBA equations for Yi^a- To see 
this let us use that Ti^i = Ti^iJ^. Then: 

f +iTr.i = T+if ^i^(J-+J-). (C.28) 
Let us consider this equation for Re (u) > 2g and shghtly above the real axis. Now let us note that 



J-+ = J- = y'T+o = y/T^ and logJ-+/J-+ = -disc(J"+) = logyi,iy2,2 - Ti,o/To,i. Since the 
T-gauge is Z4-symmetric, Tr^j^Tj^j^ = 'r2,i. On the other hand, from (3.3), (3.45) and Hirota equations 
in the magic one gets ¥2.1 — Ti^2- Plugging these relations into (C.28) we get 

t+it^i - ti,2ti^o. (C.29) 

Since the full Hirota equation for Ti_i in the magic is T^^T^^ = Ti_2Ti^o + T2,iTo,i and T2.1 7^ 
we conclude that 

To,i = . (C.30) 

Therefore we can derive equations (C.22),(C.23),(C.25) from the Z4 symmetry of the T-gauge. 
Note that this derivation is reversible and we can use equations (C.22),(C.23),(C.25) as an input to 
derive the Z4 symmetry of the T-gauge. Finally, (C. 22), (C. 23), (C. 25), (C.27) is the only input needed 
for the derivation of the formulae (D.8) from the next appendix that explicitly show that the generic 
form of the Z4 symmetry (3.23) holds. 



D Details of the Wronskian solution for the upper band 

In this appendix, the Wronskian formalism of section 4 is applied to derive the analyticity properties 
in the upper band. We also derive here the relation between various p- and q- functions and find a 
basis of these functions where all the q-functions are analytic in the upper half-plane. We also show 
that, unfortunately, there is no basis in which all q-functions of the upper band have a finite number 
of Z-cuts. 

D.l 1/2-Baxter equations 

The simplest Z4 condition Tq^i ~ implies that the columns in the corresponding determinant (4.14) 
are linearly dependent: 

q = ap++ -|-/3p-|-7p"" , (D.l) 
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where a,/?, 7 are some functions of the spectral parameter. We can easily determine these coefficients 
by computing T±i.i and T2,-i. For example'*^ 

T+1,1 = q+ A p73. = ^ ^ ^ = a+T+^ , 

^ ' P0P0 

t,_, ^ A p- ^ -q- A (^p++ + ^p) ^ -^t^l, - ^t^,. 

Therefore we get a = Tj;,i/Ti^2, 7 = t]*;i/ti,2, P = -T2,i/Ti,2. Analogously, from To.-i = a 
linear relation follows 

p = aq++ +/3q + 7q"^ (D.3) 
with the same a, /3, 7 due to the LR symmetry, and we get the following pair of "1/2-Baxter" equations: 

Ti,2q + t2,ip = t^iP+++f+ip— , 
Ti,2P + T2,iq - Tr.iq+++T+iq— . 

Excluding p (resp q) from these two equations we get a finite difference equation of the 4th order 
on p (resp q) only, which has a form of the Baxter equation for the su(4) spin chain. In what follows, 
we will make some important observations based on (D.4). 

D.2 Periodic 2-forms 

The sum and the difference of the 1 /2-Baxter equations read 

(t24 + ti,2)(q + p) = t^i(q + p)++ + t+i(q + p)— , 

(t2,i - ti,2)(q - P) = t^i(q - P)++ + t+i(q - p)"" . ^ ' ' 

These are the second order difference equations and thus the Wronskians t± = (qip)"*" A (q±p)~/Ti_i 
should be i-periodic functions (playing the same role as constants for the differential equations). To 
see this one can for instance multiply the first equation in (D.5) by (q + p)A. The l.h.s. gives zero 
whereas the r.h.s. gives precisely the condition of i-periodicity of t^. It is convenient to introduce 
instead of t± their linear combinations uj = ^^'^ ^ ~ '^^'2'^ ' 

a; = ^r^(q+ A q" +P+ Ap") , X = ^r^(q+ A p" + p+ A q") . (D.6) 
Ti,i Ti^i 

These periodic 2-forms allow to relate P(3), q(3) and q, p; 

qAw"=+P(3) , pAw"=+q(3) (D.7a) 
qAx" = -q(3) , PAX"=-P(3). (D.7b) 
The derivation is straightforward. For example, using (D.l): 

qApAp— p++ApAp— 

qAw = = = +P(3)> 

-■-1,1 -'-1,2 

qApAq qA q++ A q 

qAx = z;^ = = -q(3)- 

-■-1,1 -'-1,2 
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We remind here that, as discussed in section 4.3, the relation q0 = q(4-) = P0 = P(4) = Ti,! holds. 
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D.3 Expressions for T 

Using ttie relations from the previous section we can exclude p from the expressions for T-functions 
to get: 

Ta^2 = qi+'^'ql,""' , (D.8a) 
Tas = q'^"' A q[-"l A cj["-il , (D.8b) 
= -q[J"' A q[-;l + . (D.8c) 

Here is the derivation of the last, the most complicated formula. First, we notice that the Z4 symmetric 

Hirota equation Tq^^^Tq^^ = To,oTo,2 ~ — allows to substitute Ti_i with qo ^To,o = qo- This 
allows to rewrite the first equality in (D.6) as 

w = -^(q(2) +P(2))- (D-9) 

Explicit calculation then gives: 

T..0 = 4;' A p[-;i = -q!+;i a q'^;' + ^[-iq[t;i a j^'^^ = -qt^;' A qt-;! + ;^qg;i A p[+;i 

- - J+'^l A q[-;l + :^t[+"l = -q[+;l A q[-;l + . (D.IO) 



- q(2) "-1(2) ^ jr[+a] -"0,0 - H(2) "H(2) 

The representations (D.Sa), (D.Sb) and (D.8c) explicitly show that the expressions for Ta^s in 
terms of the Wronskian ansatz (4.14) explicitly satisfy the general Z4 property T_a,s — {~^)'''ia,s 
announced in section 4. 

D.4 Relation between p and q 

Both qi and p; satisfy the same 4*^ order finite difference equation and thus each p should be a 
linear combination of q's with i-periodic coefficients. It is easy to see from (D.7a) that in fact these 
coefficients are: 

p, = Viq, , = -uJ^kx''' , (D.ll) 
where x*^ = \t'^-'*'^Xki- For that we use that the Pfaffian of Xij is —1: 

X A X = -2-^q(2) A P(2) - -2-^To,o = -2 , (D.12) 
1,1 1,1 

where we used (4.16). Similarly, the Pfaffian of can be shown to be 1 and so 

w^^- ^ = AoJ-'r\ X^' ^ le'^'^'xa - (X^^)^-''. (D.13) 

As a consequence = 1 which implies that the inverse relation has precisely the same form q^ = Vi-''Pj 
and thus it ensures the LR symmetry. In addition one can check that Tr = 0. Hence the eigenvalues 
of y are (+1, +1, -1, -1). 
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D.5 Darboux basis 

Let us now use the freedom (4.3) in tlie choice of Q-functions to simplify the relation among p, q, p(3) 
and q(3) (D.7,D.ll). In order to preserve qg = = P0 = = Ti.i we restrict ourselves to the 
transformations with detH = 1. Under the transformation (4.3) uJij and Xij transform covariantly 
while transforms by the adjoint action: 



HujH' 



X 



(D.14) 



Note that V by itself is an i/-transforniation. In addition, since V just interchanges p and q, it leaves 
(jj and X invariant. 

Since is a skew-symmetric matrix, it is always possible to choose an iJ-transformation that 
brings ujij to the Darboux form, i.e. to the matrix with uji2 — —(^21 — ^^34 — ^^^43 1 a-nd with all 
other matrix elements being zero. The remaining freedom in iJ-transformation is the sp(4) subgroup, 
and V is its particular element. The action oiHonV generates an adjoint orbit which always passes 
through a Cartan element, hence V can be brought to the form diag(l, 1, — 1, — 1). Therefore it is 
possible to find a Q-basis in which: 
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In this canonical basis we get simple relations 
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P123 = ^3 = 

P124 = q4 = 
P234 = q2 = 
P134 = qi = 



P3 
P4 
P2 
Pi 



"qi23) 
"qi24, 

-q234, 

-qi34- 



(D.16) 



Notice that we still have a residual symmetry sp(2) x sp(2) preserving both uj and x- 

The Darboux basis is a direct analog of the 1-cut basis for the horizontal su(2)-band of the section 
4.2. It is tempting to ask whether it is possible also to find a basis of q-functions all having a finite 
number of cuts. Unfortunately, as we shell show later this appears to be not possible. 

D.6 Existence of a basis regular in a half-plane 

In this subsection we generalize the analyticity arguments, given for the right band in the main text, 
to the upper band. The situation in the upper band is a little bit more complicated and in particular 
we will see in the subsequent section that it is not possible to choose a basis with just a single cut. 
However, we will prove that there exists a basis in which q^^^-forms are regular everywhere above 
Z|j._2|_i cut and p(fe)-forms are regular everywhere below Z_|fc_2|+i cut. 



We start by writing an analog of (4.8). The idea of derivation is simple: we want to write q) 



2a-l] 



in terms of T functions with a shift ^ a and of the other q^ or pj functions with shifts independent of 
a. This is straightforward to achieve using (4.14) which we copy here for convenience: 



^a,\ — -La.-l — ^{Z) P ^ "(3) 



(D.17) 



these relations are not the same as in (4.17). 
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Indeed, we notice that expressions for T^a-2\j "^^a-i,!^ "^^a.i^^ "^a^+i^i ^^'^ independent linear 
functions of four components of q|g'^ ^' with coefficients containing only a-independent shifts. There- 
fore q(3) can be written as certain linear combinations of these functions. Explicitly: 

(p[+3] ^ p[+il ^ p[-il ^ p[-3l)qP«-il ^ _p[+i] A p[-il A p[-'lTtt'i 

+p[+3] Ap[-il ApI-^lTi'li^i 
_p[+3] Apl+ilApI-^lTl^il 
+p[+3] A p[+il A pI-^lTL^Vt'i . (D.18) 

The coefficients in front of T-functions is very easy to check by acting with ApI"', n ~ —3, —1, 1, 3 on 
both sides of the equality. 

Now we make the following observation. Let us first shift the argument in both sides of (D.18) 
by — 5i/2 so that q'g)"^' is expressed in terms of T^a-2\' "^a^-un "^^aA^^ "^a^+i^i '^i^h coefficients 
depending on p'"^', p'~^', p'~^', p'^^'. Notice that all these T-functions are shifted such that they 
are still regular on the real axis for a > 4. Thus if we find a basis where p does not have cuts 
Z-2, Z^4, Z_6, ^-8 then q(3) will be immediately regular in the whole upper half- plane. 

Note that we can reverse all shifts and interchange q(3) with *p in (D.18). From that equation 
we can say that in the basis where q does not have the cuts ^+2, .2^+4 , Z^q, Z^g, p must be regular 
everywhere in the lower half-plane. 

Using the argument from the main text we can argue that the cuts Z2, Z^, Zg, Z^ in q(3) can 
be always canceled by an appropriate H-transformation. Then in this basis we see that q(3) will be in 
fact free of all cuts above the real axis and p will have no cuts below the real axis. Note that in this 
basis P(3), as a certain combination of p^"^', p, p'^', should be regular at least starting from —i and 
below. 

To complete our proof let us make two further observations. First, due to the LR symmetry 
Ta,i = Ta,-i we can also make a replacement q — *q(3) and p — *P(3) in (D.18). Second, we notice 
that the regularity argument also works when in p(3) there are no cuts Z_4, Z^q, ■Z'-s, -Z-io instead. 
For that one should shift (D.18) by —6i/2 and use it for a > 5. Thus we conclude that in the same 
basis q and q(3) are regular in the upper half-plane and p and p(3) are regular in the lower half-plane. 
It remains to show that q(2) and P(2) are also regular in their half-planes. It is immediately obvious 
that, say, q(2) is regular starting from Zi. In fact we can prove that the cut Zi is absent and the 
function is regular everywhere above Z-i. 

Indeed, let us rewrite the Pliicker relation (4.21) in the following form: 

'iki'ikj q/cj q^i 

The l.h.s of this relation is regular everywhere above the real axis. In particular, the discontinuity at 
Z2 should be trivial, which gives for the r.h.s.: 

qfej q^r 



■^^the * is the Hodge dualization. 
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We see that o is some universal function which should be the same for any pair of indexes, since q^j 
is an antisymmetric tensor. Consider now 

T- - a+ - ^e'J'''-'a[+^lD["''^l fD 20) 

-^2,0 — 4(2) f (2) ~ Vkl ■ l^u.zu; 

Since disc (T2'o) = 0, one gets 

^i^c(qjj- jpj,; -o-e-' Pm - 013,0 ■ (U-^^j 

We see therefore that o = and hence disc (q^j) = 0, i.e. we proved that q(2) is regular above Z_i-cut. 
The same logic applies to p(2) which is thus regular below the Zi-cut. This accomplishes the proof 
about the possibility to construct a basis which is regular in a half-plane announced in the beginning 
of this subsection. 

D.7 No-go theorem: There is no q-basis with a finite number of short cuts 

Proof: Suppose there exists such a basis. Then by considerations from the previous subsection q(2) 
should be analytic in the upper half-plane up to Z_i and in the lower half-plane up to Zi, i.e. it 
should be analytic in the whole plane. The same is true for P(2)- Therefore Ta,o in (4-14) and hence 

Y„ = T^:. — 1 should not have any cuts at all which leads to an obvious contradiction with 

the basic properties of the Y-system, starting from the large L asymptotics. 

This no-go theorem shows that the Darboux basis cannot be in this respect as nice as the 1-cut 
basis of the right(left) band. It even cannot be a half-plane regular one. Indeed, in the Darboux basis 
P12 = q34 and hence, due to (D.9) IF — Tbj\i = q(2) + P(2) = qi2 + q34- IF definitely has a Z\ cut, 
hence q(2) cannot be regular at Z-y. However, the existence of a periodic form with the periodic 
structure of short Z-cuts, may potentially allow to deform the Darboux basis in a way to parameterize 
the upper band in terms of densities with the finite [—2(7, 1g\ support only. 

D.8 Basis used for computations 

In this subsection we show existence of the basis which is used for explicit computations^''. This basis 
is regular in a half-plane and posses particularly chosen relations (4.17) among p-s and q-s. 

In addition to LR and Z4 symmetries which manifest themselves in existence of respectively 2- 
forms X and lo there is also a reality property of T-functions which we will now encode into the 
properties of an i-periodic 4x4 matrix S defined by: 

qfe = iSk^Vi- (D-22) 

The latter relation should exist since complex conjugated quantities p and q satisfy the same Baxter 
equations (D.4) as p and q. It also follows from (D.4) that 

Pfc=z5fc^q,. (D.23) 

5 transforms under the i?-transformation (4.3) as S* — > HSH^^ and enjoys the following proper- 
ties: 

det5 = 1, (D.24a) 
SS = 1, (D.24b) 
SxS^ = X ■ (D.24c) 



Note that this not a Darboux basis of appendix D.5. In particular, structure of uj is complicated. 



- 52 - 



These properties follow correspondingly from (4.16), reality of T-functions, and definition (D.6) of x- 
Relations (4.17) mean a particular choice of x and S: 



Xo 
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Let us now show that this particular choice can be always done by using the ff-symmetry (4.3). First, 
we bring x to the form xo- After that we only can use the H-transformations restricted by: 



xo = HxoH' 

It is easy to construct such transformation explicitly:^^ 

H = ^/S~oyfS 

Since S transforms as S" HSH~^ we get 



S 



S ( 



= So . 



(D.26) 



(D.27) 



(D.28) 



Note that once q's are chosen to be regular in the upper half-plane (and p's are chosen to be regular 
in the lower half-plane), x and S become also regular. Hence all the transformations discussed in the 
previous paragraph are done by the regular iJ-matrices and thus do not spoil analyticity properties 
of the q-basis. 



E Further details about FiNLIE equations 
E.l Uniqueness of the physical gauge 

We are going to show that conditions (3.44) completely constrain the T-gauge, up to an overall 
constant. First, let us note that the gauge transformations that preserve the LR symmetry Tq^_s = 
Ta,s, reality Ta,s — Ta,s, and analyticity Ta.s G ^i+|a-s| are of the following form: 

Ta,s ^ ff[''+^l5f""'l5'"""'lff'""+'lT,,„ (E.l) 

where the function g is analytic in the upper half-plane above Z_i and g is the complex conjugate of 
9- 

To preserve the T2,3 — T3.2 condition g/g should be a periodic function. To preserve the T^J^q = 
Tqq condition gg should be also a periodic function. Hence g — \/gg x g/g is periodic by itself. 
Then, since g is analytic in the upper half-plane and periodic, it is analytic everywhere. Such periodic 
analytic function should not have poles, due to the condition of absence of singularities, and no zeroes, 
due to the condition of a minimal possible number of zeroes. Therefore this gauge transformation can 
be only a constant. 

Since S is an element of SP{4) (see (D.24c)) it can be written as e'^ where R is an element of the algebra sp{4) 
and R = R due to (D.24b). We define v/S = e'«/2. One can see that \/S\/S = 1 and also (D.24c) holds for v^. In the 
same way we define \/5o. 
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E.2 Reduction of Bethe equations to computable quantities 

We need to express (3.35) in terms of quantities that can be explicitly computed. We start from the 

'2,oTo,o T 

in figure 4 is trivial. Therefore (3.35) can be equivalently written as: 



fact that Yifl = x^^r — ~ t7^- ^^'^'^'^ T2,o G A2, its analytic continuation along the contour 7 defined 



T2,o(Uj) T2,o(Wj 



= -1. (E.2) 



It is allowed to replace 7 with 7''" because the branch points are of the square root type and analytic 
continuation along 7 and 7+ leads to the same results. To find explicitly the analytically continued 
function T^^ one uses Ti.i = /i+/i~7i,i and the density representation (6.3) Ti^i = 1+JCi* p. For the 
latter the continuation picks a pole in the Cauchy kernel leading to 

Ti^r-TM-P-. (E.3) 

Analytic continuation of h can be found from the discontinuity property (5.37) which gives 

{h-y^ = ^— ^^'^ = f^^' u = uj. (E.4) 

At the last step we used that 2 ^ P^^^ Bethe root. Note that F has a zero at the same 

position, therefore the whole expression is regular. 

The function p~ can be found from equation (6.4) shifted by — i/2. Again, at the Bethe root Y22 
is singular, therefore we can simplify the l.h.s of (6.4) just to (1 + Y^^)~'^ . Explicitly, we get after 
some algebra: 

'1,2 ' 1,1 ^ 1,1 

By combining (E.3), (E.4) and (E.5) we derive the final answer for Ti ^: 

= -^^y2:2i^y ^ = ur (E.6) 
'1,1 



Let us now check the following property: 

Tl^,{u,)Tl,{u,) 
T2,o(wj) 



1. (E.7) 



Applying Hirota equations on the magic sheet and using the relation between T- and T-gauges one 
can get: 

T2^0=P++P— - Til. (E.8) 

The first term on the r.h.s. can be further transformed, using the i-periodicity of J- in the mirror, the 
relation (3.41), and the Y-system equation (2.1) for Vi,!, to: 

= -F^ (11,1^2,2)+ (n,i>^2,2)" = ^^Y+^Y-^il + yi,o)(l + ri,2), u = u,. (E.9) 

Note that at the last step we used that 1 + 1/5^2, 1 = 1 at the Bethe root. 
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On the other hand, it follows from (E.6) and its complex conjugate: 

n^-^li = J''Y+,Y,:,'^^^ = T'Y,+,Y,-,il + Y,a) ■ (E.IO) 

12,2 

The difference of (E.9) and (E.IO) gives: 

Yio -^1.0-^1,0 



-^^^2!2^2,2^1,o(l + 5^1,2) - 1 - + i - -Z— — — Z - irf,i, 

J- + J-i.nJ-i.n T+^T^T 



which indeed assures (E.7). 

The property (E.7) in particular implies that if (E.2) holds then the conjugated equation 
(Ir^^(Mj))^ — ~T2,o(wj) also holds, i.e. the Bethe roots always come in the conjugated pairs. For 
the sl(2) sector all the Bethe roots are expected to be real and the eq.(E.7) follows from (E.2) up to 
a sign. We hence checked what sign should be chosen. 

The ratio of (E.2) and (E.7) leads to: 

'■'^ " = -1, (E.ll) 

nl(«,) 

^2:2^,7^'^ 

This is the equation that we are using to determine the positions of the Bethe roots. 

One can show from (5.27) and the explicit definitions for B that (E.12) is equivalent to 

(fj-V ( ^- V #[-'1 



which after substitution (E.6) gives 



(E.12) 



2.0 



(E.13) 



Then, by using explicit asymptotic expressions from appendix B it is easy to see that (E.13) reduces 
to 



Bi-)- 



' BES 



(E.14) 



in the large volume limit. This is nothing but the asymptotic Bethe Ansatz equations for the sl{2) 
sector, as it should be. 

E.3 An alternative equation on F 

The z-periodic function T = ^JTq^ plays the key role in our constructions since it is the function that 
relates the right band with the upper band. In this section we derive a simple expression for J- in 
terms of the product Yi,il2,2 which allows for a better understanding of the properties of F. For that 
we depart from (3.41) which can be considered as the Riemann-Hilbcrt equation on T . It is not hard 
to solve this equation. One of the possible solutions is''^ 



J-o{u) = exp 



^ ^ (tanh 7r(u - w) + sign(u) ) log ( Yi.i(f)5^2.2(w) 



(E.15) 



Recall that is j-periodic and hence it has an infinite set of long Z-cuts. 
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where sign(w) is added for the convergence of the integral. In general one can multiply J-'o{u) by an 
arbitrary periodic function analytic in the whole complex plane. The remaining function can be fixed 
from the knowledge that has simple zeros at the positions of Bethe roots: 

M 

J'(u) = J'o{u)Ajr Y[ sinh(7r(u - u^)) , (E.16) 
1=1 

where Ajr is a constant. 

One should be careful with the 27riZ ambiguity of log(yi_il2,2) in the definition of J-q. This 
ambiguity is however uniquely fixed by the requirement that logJ-p should not have any logarithmic 
poles at u = I ± 2g. This is only possible if the integrand is zero at v — zL2g which uniquely fixes the 
branch. 

Let us now analyze the large u behavior of J-o(u). Expansion (3.40) implies that log(Yi_iy2,2) — >■ 
2iTim± at w — > ±oo + iO. In order to determine the integers m-t we analytically continue the logarithm 
from u = 2g, where we already fixed log(Yi,il2,2) ~ 0, to ±oo + iO along the real axis. At weak 
coupling the Y-functions are known explicitly, see (B.22), and we get: 

- n ^^^^^ (E.17) 
^■^ u — Ui — 1/2 

from where we find 77i_ — m+ = M . By the continuity argument, since m± are integers, this relation 
should hold for a finite coupling^'^. For the case of symmetric roots to_ = —in^ = M /2. 

One can check that the exponential factor from J^) iu) cancels against the one from the product of 
sinh in J-{u) and we get a simple power-law asymptotics ~ , in complete agreement with (3.43). 

E.4 Fixing the normalization constants 

A freedom in the overall scale of the T-gauge was used to choose the normalization for / in (5.20). 
There is no further freedom and the integration constants A and Ajr appearing in this article are fixed 
once the normalization for / is chosen. 

Kjr appearing in (E.16) is found from the comparison of the large u asymptotics of = //■\/ >3^.o 
and (E.16) and is given by 



Kjr = CXp 



-2TiMg + E log 2g + M log 2 - / dv[~i log {Yi^i{v)Y2^2{v)) + ttM 

J2q+iO \ ' ' ^ 



(E.18) 



where the branch of log Yi^iy2_2 is the same as in (E.15). 

The comparison of (E.16) and //\/ ^o.o at large u reveals also a practical expression for the energy, 

1 

E = M - - / dvpeiv), (E.19) 



2 _ 

which is suitable for the numerical computation of E. 

A appearing in (6.12) is fixed in two steps. First, we demand that"'''^ qg = Ti i is a real function 
on the magic sheet which allows us to fix the phase of A. Indeed, by studying the equality qg = 



^^unless there is a critical point at some value of the coupling - the possibility which we ignore since it was not realized 
so far for explicit examples. 

^*q0 is regular in the upper half-plane above Zi. This is the region where q0 on mirror and magic sheets are identified. 
In contrary, mirror Ti i and magic Ti i are identified inside their analyticity strip ^i. Hence q0 and Ti i are equal in 
magic but not in mirror. 
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%f^^ f U^U at large u we see that is a real function, at least at large u, if A^e**'''= is a real 
function at large u. Hence: 



which should be valid for —2g < u < 2g. This relation is an upper band analog of (5.37) and it is 
derived similarly. 

Note that (E.21) fixes UU on Zq essentially from the knowledge of p2. Hence the nontrivial 
information about U is contained in the complementary region: {—oo,—2g] U [2g,cxi). 
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